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INTRODUCTION AND SUMMARY 

Interstellar dust is an important component of the interstellar 
medium. It dominates the heating and cooling of clouds through ener- 
getic photoelectrons and gas-grain collisions (see the chapter by Black 
in this volume). It is also the dominant source of opacity and deter- 
mines, therefore, the spectrum of dust-enshrouded objects. Unquestion- 
ably, molecular hydrogen, the most abundant gas-phase species in molecu- 
lar clouds, is formed on grain surfaces. Other gas-phase molecules 
might result from reactions on grain surfaces as well. Grains can also 
influence the gas-phase composition of molecular clouds indirectly, 
because they may lock up some of the elements. The metals (e.g., iron) 
are particularly important in this respect since they tend to carry the 
charge inside a dense cloud and, thus, regulate the ion-molecule chemis- 
try (see the chapter by Prasad in this volume). Such depletions will, 
of course, also influence the cooling balance. 

Yet, despite the importance of the dust and despite 50 years of 
active research in this field, many of the major questions either remain 
unanswered or the answers that have been proposed are highly controver- 
sial. These include the two basic ones: What is the composition of the 

interstellar dust? and Where and how is it formed? These questions are 
of course interrelated and an answer to one suggests an answer to the 
other. 

The dust components proposed to be present in the interstellar 
medium can be divided into two classes depending on their formation 
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history. First, there is stardust, such as silicates, graphite and 
amorphous carbon, which are made in the oxygen-rich or carbon-rich 
outflow from late-type giants and planetary nebulae and, possibly, novae 
and supernovae. Dust formation has been theoretically predicted to 
occur in such outflows (Hoyle and Wickramasinghe 1962; Kameijo 1963). 
This was subsequently confirmed, mainly by infrared observations, in the 
early seventies. There is now abundant evidence for the formation of 
dust in these high-density, high-temperature environments (cf. the 
reviews by Merrill 1977 and Aitken 1981). Second, there are dust compo- 
nents that are formed in the interstellar medium itself. This includes 
"icy" grain mantles, consisting of simple molecules (e.g., H2O, NH^, 
CH^OH, and CO) inside dense clouds, as well as an organic refractory 
dust component, consisting of more complex molecules, in the diffuse 
interstellar medium. The presence of both of these dust components in 
the interstellar medium has been revealed by infrared spectroscopy of 
highly obscured objects. The icy grain mantles, which form by accretion 
and reaction of gas-phase species on preexisting cores, seem to be 
ubiquitous inside molecular clouds (Merrill et al. 1976; Joyce and Simon 
1982; Willner et al. 1982). In contrast, the organic refractory compo- 
nent has, at present, only been detected toward the galactic center 
(Willner et al. 1979; Allen and Wickramasinghe 1981; Butchart et al. 
1985). Presumably, there is a chemical, evolutionary link between these 
two dust components. That is, the organic refractory dust component is 
formed by energetic processing (UV photolysis or cosmic ray bombardment) 
of the icy grain mantles (Greenberg 1982; Moore and Donn 1982; Strazulla 
et al. 1983) . 

Detailed models of interstellar extinction have been constructed 
based on stardust alone (e.g., Mathis et al. 1977; Draine and Lee 1984), 
as well as on a combination of stardust and interstellar medium dust 
(Greenberg and Hong 1974; Greenberg and Chlewicki 1984). Not too much 
credence should be given to such models since the interstellar extinc- 
tion curve is quite insensitive to the exact composition of the dust. 
This is because other parameters (e.g., size distribution, shape, and 
surface roughness) will influence the extinction curve of an ensemble of 
particles as well. This non-uniqueness is compounded by the variations 
in the optical properties, even within one class of materials such as 
silicates, and the resulting uncertainty in the optical properties of 
interstellar grain materials. At this moment, the most diagnostic 
information on the dust composition results, therefore, from features in 
the extinction curve, notably in the infrared, rather than from the 
overall shape of the extinction curve. 

In §1, a compilation of the different dust components observed to 
be present in the interstellar medium and their relative importance is 
presented. These include the stardust components amorphous silicates, 
amorphous carbon, polycyclic aromatic hydrocarbons (PAHs), and graphite, 
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as well as the interstellar medium components , which are the organic 
refractory grains and icy grain mantles. These data are summarized in 
Table I. It is emphasized that the uncertainties in these observations 
preclude a definite choice between the different proposed dust models, 
although it seems that stardust may dominate the dust composition in the 
solar neighborhood. 

In §2, the difference between amorphous and crystalline materials 
is examined. Some attention is given to the structure of truly amor- 
phous materials (e.g., random network materials) such as glasses (e.g., 
silicates) and ices. The microcrystalline structure of soot particles 
is also discussed with an emphasis on the structural differences and 
similarities between graphite sheets, amorphous carbon particles, and 
PAH molecules. This section is concluded with a discussion of the 
structure of amorphous polymer solids, appropriate for the organic 
refractory dust component. 

In §3, the optical properties of amorphous and crystalline mate- 
rials are contrasted. The emphasis is on silicates (amorphous), because 
the properties of carbon grains are to some extent already covered in 
the chapter by Allamandola, Tielens, and Barker in this volume. The 
optical properties of materials divide naturally into three regions 
which are discussed separately. The first is the optical and UV region, 
which is dominated by the electronic properties of the material. Sec- 
tion 3.1 focuses on understanding the observed, enhanced absorption at 
about 2 urn of silicates formed in red giant envelopes. This absorption, 
which lies within the band gap of silicates, is due to localized elec- 
tronic states, probably caused by the presence of impurities such as FeO 
in the silicate. The second is the mid-infrared region, which for 
amorphous materials is dominated by the fundamental vibrational modes of 
the nearest-neighbor atoms. Special emphasis is given in §3.2 to the 
large width of the 10 um Si-0 stretching vibration in amorphous sili- 
cates. It is suggested that the decreased width of the 10 um feature in 
absorption as compared to- emission is due to the low temperature of the 
dust in the interstellar medium. The third is the far-infrared region, 
which shows the most distinctive differences between amorphous and 
crystalline materials. In §3.3 the conditions for a X" 1 far-infrared 
absorption efficiency are examined. From interstellar far-infrared 
observations it is concluded that the volume of the interstellar dust is 
dominated by layered materials, such as amorphous carbon or layer- 
lattice silicates. 

The physical principles of grain surface chemistry are discussed in 
§4. This includes the surface structure of materials (§4.1), the 
adsorption energy and residence time of species on a grain surface 
(§4.2), the sticking probability (§4.3), the mobility of species on a 
grain surface (§4.4), the reactions that are of importance for the 
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formation of icy grain mantles (§4.5), and the dissipation of the reac- 
tion heat (§4.6). 

It is generally agreed that H 2 is formed on grain surfaces in tht 
interstellar medium. Because of its low mass, a physically adsorbed 
H atom can tunnel through the activation barrier for reaction with a 
chemisorbed H atom and form H 2 . This process is examined in some detail 
in §5. 

In §6 the depletion of elements due to their accretion on grain 
surfaces in the diffuse interstellar medium i3 examined. It is sug- 
gested that the high depletion of elements such as Fe and Ca is due to 
chemisorption on grain surfaces. Other atoms (e.g., 0 and N) will also 
be chemisorbed, initially, but reaction with H will form saturated 
hydrides which are easily photodesorbed. Their depletion will then be 
fairly low. 

The composition of interstellar grain mantles is discussed in §7. 

It results from a complex interplay of gas-phase and grain-surface 
reactions. Broadly speaking, three different regimes can be discerned, 
depending on the conditions in the gas phase. First, a reducing atmo- 
sphere in which most of the accreting gas-phase species are saturated 
through reactions with atomic hydrogen on the surface. Second, an 
oxidizing atmosphere in which surface reactions with atomic oxygen 
dominate. Third, an inert atmosphere in which the accreting gas-phase 
species essentially do not react with each other. Observations show 
that the reducing conditions dominate during the formation of inter- 
stellar icy grain mantles. 

In §8 the formation of the refractory dust component that is 
observed in the diffuse interstellar medium is discussed. Essentially, 
the problem is to get from the simple, largely hydrogenated molecules in 
icy grain mantles to more complex organic refractory material, which is 
able to survive better in the harsh environment of the diffuse inter- 
stellar medium. Some form of energetic processing of the grain mantle 
is required for this transformation. The most likely sequence of events 
in UV photolysis producing radicals, followed by the diffusion and 
reaction of these radicals to produce complex molecules. In order to 
get a substantial yield of organic refractory dust material, it is 
concluded that there has to be a significant internal UV field in molec- 
ular clouds and that transient heating mechanisms such as cosmic ray 
heating, which promote diffusion, are important. 

Finally, in §9 the contribution of grains to the gas-phase composi- 
tion of molecular clouds is examined. It is emphasized that the criti- 
cal temperature for grain-mantle formation is only about 30 K and that 
above this temperature grain-mantle formation is a slow process. Ejec- 
tion mechanisms from grain surfaces are discussed in §9.1. Transient 
heating of photolyzed grain mantles by cosmic rays or grain-grain colli- 
sions up to about 30 K promotes diffusion and reaction of stored radi- 
cals. The liberated reaction heat will explosively evaporate an 



appreciable fraction of the icy grain mantle. During this process more 
organic refractory dust material can be produced. The signatures of 
grain chemistry are discussed in §9.2. It is pointed out that grain 
chemistry will mainly lead to simple hydrogenated molecules (e.g., 

H 2 0, NHj, and CHgOH), large deuteration effects (e.g., 

a 1:0.1:0.01:0.001), saturated metal hydroxides (e.g., 
NaOH), and ortho/para ratios characteristic for the dust temperature. 


I. COMPOSITION OF DUST IN DIFFUSE INTERSTELLAR MEDIUM 

Many different dust components have been proposed to be present in the 
diffuse interstellar medium. These include silicates, amorphous carbon, 
PAHs, graphite, organic refractories, SiC, metallic oxides and many more 
(cf. Mathis 1986 and references therein). Most models rely on some 
combination of these to explain the interstellar extinction curve. 

Direct information on the dust components present in the diffuse inter- 
stellar medium comes, actually, from structure in the extinction curve; 
for example, the 2200 k bump, the 10 and 20 um features; several infra- 
red absorption features in the spectrum of the galactic center; a 7.6 um 
absorption feature in some protostars; and the 3.3, 6.2, 7.7, and 

II. 3 um emission features. A successful model should, of course, also 
explain the wavelength dependence of the observed extinction curve and, 
in particular, the amount of visual extinction per hydrogen atom 
A V /N H . These constraints can be used to identify the dust components 
present in the diffuse interstellar medium and to determine their rela- 
tive importance. In the remainder of this section we present such an 
analysis of these constraints. It should be emphasized, however, that 
not all researchers in the field of interstellar dust will agree with 
our interpretation of the data, in particular with the details. The 
reader is referred to the review by Mathis (1986) for a critical analy- 
sis of some of the different dust models proposed. 

The spectroscopic constraints listed above can be used to identify 
the different dust components present and to determine their relative 
importance. In line with most other analyses, we attribute the ubiqui- 
tous 2200 k bump and the 10 and 20 inn features to graphite and silicate, 
respectively. Evidence for the presence of amorphous carbon, organic 
refractory material, and PAHs is provided by a 7.6 um absorption fea- 
ture, several other absorption features, and the IR emission features, 
as detailed below. In Table I we summarize the estimates of the mass 
and volume per H atom, Mj and Vj, of each dust component identified in 
the interstellar medium. These J have been determined from the observa- 
tions using 
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Table 1: The Composition of the Dust in the Diffuse Interstellar Medium 
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cm - /H atom. 


and 


Vj = y » 3 , (?) 

where f^ is the fraction of element i locked up in dust component 
J; is the solar abundance of element i; m^ is the atomic mass of 

element i; o s is the specific weight of the dust component; and where 
the summation is over all elements in a dust component. The fraction of 
an element locked up in a particular dust component (relative to solar 
abundances) is estimated from the observed strength of the absorption 
(or emission) features per hydrogen atom. The details of this analysis 
are given below. 

1.1 Silicates 

The 10 and 20 um features are seen in many objects in absorption as well 
as emission, including oxygen-rich M giants (presumably the birth 
sites), HII regions, protostars, and background stars (Merrill 1977; 
Aitken 1981). Generally, these features are attributed to the Si-0 
stretching and bending vibrations in silicate materials. The width of 
the observed features (e.g., much broader than in terrestrial minerals) 
and the absence of substructure suggests that the silicates have an 
amorphous structure or are hydrated (Zaikowski et al. 1975; Day and Donn 
1978; Day 1979). 

The fraction of the silicon locked up in silicate dust grains in 
the diffuse interstellar medium can be estimated from the observed 
strength of the 10 um silicate absorption feature per magnitude of 
visual extinction. For eight bright, nearby (<3 kpc) WC Wolf-Rayets 
and supergiants, the ratio t (10 ym)/A v is determined to be 
18.5 + 1.5 mag -1 (Roche and Aitken 1984). Assuming a standard dust-to- 
gas ratio, N H /A y , appropriate for the diffuse interstellar medium of 
1.9 * 10 21 mag”* cm" 2 (Bohlin, Savage, and Drake 1978) and using the 
observed width of the 10 um feature of 220 cm” 1 , this translates into an 
integrated strength, W (10 um), of 6.2 x lO”^ 1 cm per hydrogen atom. 
Because of the disordered nature of amorphous or hydrated silicates, it 
is expected that all modes will be infrared-active and, thus, that the 
integrated strength of the Si-0 stretching vibration will not vary much 
from one silicate material to the other (see §3.2); this is borne out by 
laboratory studies. Laboratory measurements of the integrated strength 
of the 10 um feature in amorphous or hydrated magnesium silicates typi- 
cally yield 1.2 x 10” 1 ® cm per Si atom (Penman 1976; Day 1979; Mooney 
and Knacke 1986). Note that the measured peak strength and width vary 
by a factor of 3, but the integrated strengths agree to within 30%. 

Using this value requires about 150% of the solar Si abundance to be 
locked up in silicates. Such an elemental enrichment in the present 
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interstellar medium is perhaps not impossible. However, larger 10 um 
integrated absorption strengths have been reported for crystalline 
silicates (2.8 * ICT 1 ^ cm per Si atom; Zaikowski, Knacke, and Porco 
1975). Using this value requires that at least 63% of the silicon be 
locked up in silicate dust grains in the interstellar medium. Quite 
likely, the actual fraction of Si locked up in silicates is about 100%, 
and that is what we have assumed in Table I. Based on cosmic abundances 
and condensation sequences (in thermodynamic equilibrium; thus, perhaps 
not completely applicable) a magnesium or iron-magnesium silicate is 
expected to condense out in the outflow from late-type giants (Hackwell 
1971). Some Na, Ca, and Al might also be contained in interstellar 
silicates. A composition of (Mg,Fe)Si04 has been assumed in the 
calculation of the silicate mass per hydrogen atom. The specific 
density of terrestrial silicate materials ranges from about 2.5 to 
3.5 g/cm^. Interstellar silicate materials are probably amorphous and 
their specific density is on the low side of this range. In the calcu- 
lation of the silicate volume per hydrogen atom (Table I) a specific 
density of 2.5 g/cm^ has been assumed. Note that both these assumptions 
tend to maximize the contribution of silicates to the interstellar dust. 

1.2 Amorphous Carbon and Polycyclic Aromatic Hydrocarbons 

A broad and shallow 7.6 um absorption feature has been observed in the 
spectrum of NGC 7538-IRS9 and W33A, two protostars with the largest 
known 10 um optical depth (Tielens et al. 1986a). This band may be due 
to the CC stretching vibration in amorphous carbon dust particles, 
although other interpretations (in terms of molecules in icy grain 
mantles) are also possible. The hydrogen column density or visual 
extinction toward these sources is not known, but from a comparison with 
the 10 um optical depth and assuming that all of the Si is in the form 
of silicate grains, the fraction of the carbon locked up in amorphous 
carbon grains is calculated to be 5-10% (Tielens et al. 1986a). If the 
silicon depletion in silicate grains is actually less, then this number 
has to be revised downward accordingly. The resulting mass and volume 
estimates (assuming a specific weight of 2 g/cm^) for amorphous carbon 
are given in Table I. 

Emission features at 3.3, 6.2, 7.7, and 11.3 uffl have been observed 
in a variety of sources, including planetary nebulae, HII regions, 
reflection nebulae, and galactic nuclei (Russell, Soifer, and Willner 
1977; Aitken and Roche 1984; Sellgren, Werner, and Dinerstein 1983; 

Cohen et al. 1986). These features have been attributed to vibration 
involving polycyclic aromatic hydrocarbons in various forms (Duley and 
Williams 1981; Leger and Puget 1984; Allamandola, Tielens, and Barker 
1985). From an analysts of the emission mechanism, the fraction of 
carbon locked up in these molecules is estimated to be about 1%. For 
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completeness, this component is included in Table I. It is, however, 
likely that this component is just the molecular domain of the size 
distribution of interstellar amorphous carbon grains (see §3.5). 

1.3 Organic Refractory Material 

Absorption features due to an organic refractory dust component have 
been detected along the line of sight toward the galactic center 
(Willner et al. 1979; Allen and Wickramasinghe 1981; Butchart et al. 
1985; Tielens et al. 1986b). The dust along the line of sight toward 
the galactic center is presumably located in the diffuse interstellar 
medium (Roche and Aitken 1985). From an analysis of these spectra in 
terms of functional groups of simple molecules, the fraction of the 
carbon and oxygen locked up in this dust component has been determined 
to be 24J and 6J, respectively (Tielens et al. 1986b). The H column 
density has been determined from the estimated visual extrinction, 

A v = 30 mag, assuming a standard dust-to-gas ratio. The specific den- 
sity of this dust component has been assumed to be 1 g/cm^. 

It should be emphasized that infrared spectra generally only yield 
the characteristic molecular subgroups responsible for the absorption 
and not the specific molecule. The intrinsic absorption strength of 
these molecular subgroups, and thus the calculated column density, is 
therefore somewhat uncertain (by a factor of 2). Furthermore, it should 
be noted that the silicate abundance in the galactic center is larger 
than in the local interstellar medium (Roche and Aitken 1985); this may 
also pertain to the organic refractory component. In fact, the 3 um 
spectrum of the supergiant VI Cygni No. 12 (-3 kpc) does not show the 

3.4 um absorption feature (r (3.4 um) < 0.03; A y s 10 mag; Gillett 

et al. 1975a). Assuming that the composition of the organic refractory 
dust in the local interstellar medium is similar to that along the line 
of sight to the galactic center, we arrive at an upper limit on the 
fraction of the carbon locked up in the dust component of about 10£. 
Clearly, the volume and mass values given in Table I should only be 
considered as order-of-magnitude estimates. 

1.4 Graphite 

The strong 2200 A bump, which is a ubiquitous component of the inter- 
stellar extinction curve in the diffuse interstellar medium (cf. the 
review of Savage and Mathis 1979 and references therein), is generally 
attributed to small (-200 A) spherical graphite grains (Gilra 1972). 
About 25 % of the available carbon is required to explain the strength of 
this feature (Tielens and de Jong 1979). Note that larger carbon grains 
show no detectable absorption features, and their abundance can pres- 
ently only be guessed at. A specific density of 2.2 g/cm^ has been 


9 



assumed for graphite. There are some problems in attributing the 2200 A 
bump to graphite particles. The graphite transition is intrinsically 
very strong and, because of the surface polarization charge, the shape 
and peak position are very sensitive to shape and size (an excellent 
review of the effects of the surface polarization charge on the spectrum 
of small particles is given by Bohren and Huffman 1984). The observa- 
tions imply 200 A spherical graphite grains. Moreover, there may be 
some problems with graphite in explaining the observed variations in the 
ratio of the FUV extinction (X < 1 600) to the 2200 A bump (Greenberg and 
Chlewicki 1984). For these reasons, the 2200 A bump has sometimes been 
attributed to amorphous carbon instead of graphite grains (cf. Mathis 
1986). (The structural difference between amorphous carbon and graphite 
will be discussed in §3.5). Indeed, the intrinsic strength of the 
amorphous carbon transition is much less (Borghesi et al. 1983), and 
surface charge effects are, therefore, unimportant. The drawback of 
amorphous carbon is, however, that more material is required to explain 
the observed strength of the 2200 A feature per hydrogen atom. This 
does not only conflict with the amount of amorphous carbon determined 
from the infrared, but actually requires more than the cosmic abundance 
of carbon. 

1.5 Total Extinction 

A rough estimate of the total volume of dust material required to 
explain the observed extinction can be made from a Kramers-Kronig analy- 
sis of the extinction curve (Purcell 1969). An extensive discussion of 
this analysis is given in Spitzer (1978). For reasonable values of the 
optical properties of the dust, the value in Table I results. This is 
only somewhat larger than the sum of the identified components. A more 
detailed analysis of the extinction curve is, however, required to 
determine whether the identified components can indeed explain the full 
extinction curve. 

1.6 Interstellar Grain Mantles 

For completeness we now mention interstellar grain mantles. Inside 
dense molecular clouds, interstellar grains can acquire an icy grain 
mantle, consisting of such simple molecules as H 2 O, CHgOH, CO, and NHg 
(Merrill et al. 1976; Hagen et al. 1980; Knacke et al. 1982; Lacy et al. 
1984; Tielens et al. 1984, 1986a). The fraction of the carbon and 
oxygen locked up in these grain mantles ranges typically from 1-40% with 
a corresponding mass of 0.2-10 x 10" 2 ^ g per H atom and corresponding 
volume of 0.2-10 x 10 -2 ^ cm^ per H atom. Thus, the growth of these icy 
grain mantles in molecular clouds can actually double the total dust 
mass or volume. However, these icy grain mantles have never been 
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observed outside of molecular clouds. Presumably, these volatile mate- 
rials are easily destroyed by interstellar shocks and UV photodesorption 
and their abundance in the diffuse interstellar medium is, therefore, 
low (Barlow 1978; Draine and Salpeter 1979); as a result, they are not 
included in Table I . 

1.7 Discussion 

In the preceding paragraphs different dust components in the diffuse 
interstellar medium have been identified (cf. Table I). These are 
silicates, amorphous carbon, PAHs, organic refractory materials, and 
graphite. From Table I, one may get the impression that the organic 
refractory dust component and silicates dominate the composition of the 
interstellar dust rather than graphite. Greenberg and his co-workers 
(Greenberg and Hong 1974; Greenberg and Chlewicki 1984) have proposed 
such a model for the interstellar dust. It should be borne in mind, 
however, that the estimate of the total contribution of the organic 
refractory component is actually quite uncertain; it may be much less 
than estimated, particularly in the solar neighborhood. If that is the 
case, the bulk of the dust has to be in a form without detectable IR or 
UV absorption features. Large graphite particles fulfill these condi- 
tions. The electronic intraban'd transitions in graphite screen the 
(weak) lattice vibrations of the carbon atoms and make them presently 
undetectable (Draine 1984). Furthermore, if their size is large enough 
( >0. 1 urn) large graphite particles will have an extinction cross section 
in the UV equal to the geometric cross section and, therefore, will not 
show a detectable UV feature. Thus, if the volume of the organic 
refractory component is much less than estimated in Table I, then large 
graphite particles without detectable absorption features (a theorist's 
dream!) may dominate the composition of the interstellar dust. Models 
based on this premise have been developed by Mathis and his co-workers 
(Mathis, Rumpl, and Nordsieck 1977; Mathis and Wallenhorst 1981) and by 
Draine and Lee (1984). Choosing between these two possible models is 
presently not possible. As emphasized before, the most uncertain compo- 
nent, besides graphite, is the organic refractory material. Determining 
its composition and, hence, its volume will require studies of the 
infrared absorption features along many different lines of sight, in 
particular in the local diffuse interstellar medium. Nevertheless, it 
should be emphasized that all the materials identified in Table I are 
important components of the interstellar dust and contribute substan- 
tially to the interstellar extinction curve in at least some part of the 
spectrum. In our opinion, it is quite likely that the composition of 
the interstellar dust lies somewhere in between the two extremes repre- 
sented by the models of Greenberg and Mathis. 
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2. THE DISORDERED STRUCTURE OF INTERSTELLAR DUST 


Interstellar grains are generally thought to be highly disordered. This 
assessment is based on two observations. First, the width of the 10 um 
silicate feature is much larger than that in crystalline silicate mate- 
rials. This broadening can result from the disorder present in amor- 
phous or hydrated silicates (Zaikowski et al. 1975; Day and Donn 1978; 
Day 1979). Second, the far-infrared cross section of the interstellar 
dust seems to fall off much slower with wavelength than expected for 
crystalline solids. Laboratory studies suggest that this is also a 
characteristic of some amorphous materials. Theoretically, disordered 
structures are expected because grain condensation generally occurs 
under nonequilibrium conditions. In this section we will examine the 
structure of interstellar grains with emphasis on amorphous silicates. 

A detailed discussion on the structure of amorphous materials and their 
properties can be found in the Ehrenreich report (1972) and in the 
monograph by Zallen (1983). 

2.1 Structure of Amorphous Solids 

Amorphous solids, sometimes referred to as glasses, are characterized by 
complete chemical coordination (i.e., few dangling bonds), as are crys- 
talline solids. Similarly, the nearest-neighbor distances (i.e., bond 
lengths) are nearly constant. There is, however, a considerable spread 
in bond angles. Because of this latter point, amorphous solids do not 
have long-range order, in contrast to crystalline solids. This struc- 
tural difference between amorphous and crystalline solids is illustrated 
schematically in Figure 1. In the crystal, the nearest-neighbor equi- 
librium distances and the bond angles are all equal. Consequently, it 
possesses long-range order and translational symmetry. The atomic 
arrangement in an amorphous solid is a continuous random network that 
lacks symmetry and periodicity (Zachariasen 1932). Specifically, the 
nearest-neighbor distances are only almost equal. Moreover, there is 
considerable spread in the bond angles. Clearly, there is a consider- 
able amount of local order but a complete absence of long-range order in 
the amorphous state (cf. Fig. 1). The presence of the short-range order 
is, of course, Just a consequence of the chemical bonding of the solid. 

Although not strictly applicable to interstellar solids, it may be 
advantageous from a conceptual point of view to study a schematic 
volume-temperature phase diagram (Fig. 2). When the temperature is 
lowered, a gas will condense to the liquid phase at the boiling tempera- 
ture, T b . Essentially, this is a transition from a regime where the 
intermolecular forces are unimportant to one in which they dominate. 
Although the free-volume (i.e., the volume available to each molecule) 
is much less in the liquid phase than in the gas phase, it is still 
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Figure 1. Atomic arrangement in a crystalline and an amorphous solid. 
The dots indicate the equilibrium positions about which the atoms 
vibrate, and the lines indicate the chemical bonds. 


sufficiently large to allow translational motion to the molecules. A 
molecule in the liquid phase is continuously making and breaking bonds 
with its neighbors while it is diffusing through the liquid. This is, 
of course, what gives a liquid its macroscopic fluid character. Thus, 
because of the chemical bonding of the molecules, a liquid will have 
considerable short-range order. But, because of the continuous rear- 
rangement of the molecules, the bond angles and bond lengths will show a 
broad distribution around the average value. Consequently, there is no 
long-range order in a liquid. 

When a liquid is very slowly cooled to the freezing temperature, it 
will transform to the crystal phase (route 1 in Fig. 2). That is, 
crystallization centers will form and will grow at the crystal-liquid 
interfaces. Like the gas-liquid transition, this process is accompanied 
by an abrupt change in the free-volume available to the molecules.. In 
the crystal, the molecules are arranged on a translationally periodic 
lattice, which exhibits long-range order. It is important to realize 
that the crystallization process is essentially a kinetic one. That is, 
during the nucleation and growth the molecules have to "find" each other 
and arrange themselves on the crystal lattice in the energetically most 
favorable orientation. This takes time. Thus, if the liquid is rapidly 
cooled below its freezing temperature nucleation cannot take place and 
the liquid phase will persist (i.e., it will form a supercooled liquid). 
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Figure 2. Volume-temperature phase diagram. There are two ways for the 
condensation of a solid to proceed: route 1 is the path to the crystal- 
line state; route 2 is the rapid-quench path to the amorphous state. 


When a fluid is cooled below its freezing temperature while nuclea- 
tion is prevented (i.e., rapid cooling; route 2 in Fig. 2), the free- 
volume available to its molecules will shrink further and further 
because of the thermal contraction (i.e., the anharmonicity of the 
intermolecular potential). It will finally reach a point where the 
available free-volume does not allow any further translationary motions 
of the molecules and the fluid has been transformed to an amorphous 
solid. The molecules are "stuck" together in a rather random way. The 
molecules will now only show oscillatory motion around a well-defined 
equilibrium position. The temperature at which this transformation 
occurs, T a , depends slightly on the cooling history of the solid, 
because the rate of kinetic processes depends on the temperature. In 
contrast to the gas-liquid and liquid-crystal transformations, this 
transition is characterized by a smooth volume change. Because of the 
decreased free-volume of the molecules the distribution of the bond 
lengths in the amorphous state will be much less broad than in the 
liquid state, resembling more the crystalline state. The distribution 
of bond angles however is very broad, similar to that in the liquid 
state. Consequently, although there is short-range order, there is no 
long-range order . 

It is thus important to realize that any liquid or gas when cooled 
rapidly enough (i.e., nucleation prevented) will form an amorphous 
solid. This amorphous state is generally metastable with respect to the 
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crystalline state, but the latter is kinetically inaccessible. The 
detailed structure of a particular amorphous material will depend on the 
details of its formation. In particular, if the cooling is very rapid, 
the amorphous solid may retain the local structure of the liquid (e.g., 
the presence of highly strained bond angles). In contrast, when the 
cooling is somewhat slower, the amorphous solid may be able to rearrange 
itself somewhat and attain a thermodynamically slightly more favorable 
structure in the amorphous state. It should be emphasized that there is 
not just one amorphous state but many, each with a slightly different 
energy. In fact, there is an almost continuous distribution of amor- 
phous states with energies up to kT a , where T a is the temperature for 
the liquid-amorphous solid transition. Thus, although amorphous solids 
may approach an "ideal" metastable noncrystalline state, they can be 
quenched into many other metastable states. The structure and proper- 
ties of an amorphous solid depend, therefore, also on the preparation 
conditions. 

Such a distribution of amorphous states will be even more prominent 
for the condensation conditions prevailing in interstellar space. In 
that case, condensation will take place directly from the vapor state to 
the solid state. Often the temperature of the system will be much less 
than the transition temperature, T a , from the amorphous to the crystal- 
line state (note that in this case, this transition does not go through 
the liquid phase and the temperature T a differs from that defined 
above). Accreting molecules will, thus, have insufficient energy to 
arrange themselves favorably and the resulting structure will be far 
from that reached by cooling down from the energetically more favorable 
structure of an amorphous solid deposited at a temperature around T a . 
Upon warm-up, an amorphous structure rapidly deposited at a low tempera- 
ture can rearrange itself somewhat to an energetically more favorable 
amorphous structure. This annealing process will continue until the 
temperature at which it transforms into the crystalline phase is 
reached. This annealing process for the interstellar dust has been 
studied experimentally in considerable detail for amorphous ice 
(Hagen, Tielens, and Greenberg 1981) and in somewhat less detail for 
amorphous silicates (Day and Donn 1978; Nuth and Donn 1982). 

2.2 Impurities and Defects 

Up to now we have only discussed disorder caused by condensation below 
the amorphous-crystalline transition temperature. Disorder can, how- 
ever, also be due to the presence of impurities and defects in the 
crystal. Impurity atoms or molecules will have slightly different 
sizes, polarizabilities, bond strengths, and angles than the host 
species. Their presence will necessitate some degree of local rear- 
rangement of the crystal to accommodate them. This is, of course, 
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particularly true for an impurity ion with a valence different from the 
host in an ionic crystal. Impurities will, thus, distort the surround- 
ing host species (slightly) and thereby destroy the long-range order of 
the crystal. For icy grain mantles formes in interstellar molecular 
clouds this type of disorder is very important because these mantles 
consist of about equal parts of H 2 O and CH^OH (or another alcohol) and 
traces of other molecules such as CO and NH^ (Helens et al. 1986 a). Of 
course, such icy grain mantles are formed at temperatures (-10 K) much 
below T a and would, therefore, be highly disordered even if they had a 
homogeneous composition. 

Crystal lattices will also exhibit defects, owing to the presence 
of vacancies in lattice positions (Schottky defects) or the promotion of 
molecules or ions from substitutional to interstitial sites (Frenkel 
defects). The formation of these defects in a crystalling solid has an 
appreciable activation barrier. In interstellar space they might result 
from cosmic ray bombardment. A concentration of point-defects may give 
rise to dislocations (screw or edge dislocations). Dislocations are 
generally a product of the nucleation and growth processes of crystal- 
line solids. In fact, a crystalline solid generally grows at these 
dislocations. For an amorphous solid, formed under supercooled condi- 
tions, dislocations are, however, of little importance. 

2.3 Amorphous Silicates 

Silicate glasses (e.g., fused quartz) are the prototype of (binary) 
amorphous solids, and their structure is well represented by a continu- 
ous random network in three dimensions (the three-dimensional analog of 
Fig. 1; Zachariasen 1932). In fused quartz, each silicon atom is pre- 
dominantly covalently bonded to four nearest-neighbor oxygen atoms, and 
each oxygen atom is bonded to two nearest-neighbor silicon atoms. Each 
of the SiO^ units consists of an almost regular tetrahedron, and most of 
the structural variations are due to variations in the Si-O-Si bond 
angles (Bell and Dean 1972). 

The presence of electropositive elements such as Mg changes this 
picture somewhat. There are now two kinds of oxygen atoms: bridging 

oxygens, which link two tetrahedrons together, and non-bridging oxygens, 
which belong only to one tetrahedron. These non-bridging oxygens are 
bonded to the electropositive element by charge ransfer and are topolog- 
ically "dead-ended." For a small number of cations (positive ions), the 
resulting structure is that of a polymer with long chains cross-linked 
at intervals and cations in holes in between the oxygen tetrahedrons. 
This is illustrated in Figure 3a. When the number of bridging oxygen 
atoms becomes less than two per tetrahedron, spatial coherence by cross- 
linking is not possible. The silicate structure consists now of iso- 
lated chains (cf. Fig. 3b). If the number of cations present increases 
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even further, full charge separation may take place and each tetrahedron 
is only bonded to cations. In any case, however, the disorder in sili- 
cate glasses is due to the variable orientation of the SiO^ tetrahedrons 
with respect to each other. Additional disorder can be caused by the 
presence of different cations in the structure. 



O BRIDGING OXYGEN 
9 NON-BRIDGING OXYGEN 
@ CATION 


Figure 3. Two-dimensional representation of the atomic arrangement in 
silicate glasses. Each silicon atom is in the center of an oxygen 
tetrahedron. Some oxygen atoms bridge two tetrahedrons whereas others 
are bonded to cations. Depending on the number of cations present the 
SiOty chains will be cross-linked (a) or not (b). 


2.4 Amorphous Ice 

The structure of amorphous ice is very similar to that of silicates 
discussed above (Madden et al. 1978). Its infrared spectrum behaves 
very similar to that of silicates except for effects owing to hydrogen 
bonding. This has been discussed extensively in an astrophysical con- 
text by Hagen, Tielens, and Greenberg ( 1983a, b). Here we will merely 
examine the importance of the thermal history for the structure of 
interstellar amorphous ice. Consider the condensation of H 2 O molecules 
on preexisting silicate cores in the outflow from oxygen-rich late-type 
giants. For the typical vapor pressures of H 2 O in these outflows, the 
condensation temperature of ^0 is about 90 K (Nakagawa 1980; that is, 
at this temperature the accretion rate equals the evaporation rate). 

The time-scale for the amorphous-crystalline transformation is given by 
a Boltzman expression: 
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( 3 ) 


t = v ” 1 exp(E /kT) , 
a 

where E a is the activation energy for the transformation and v is a 
characteristic frequency. Because the disorder in amorphous HjO ice is 
mainly a result of a broader distribution of bond angles (Madden et al. 
1978), it is appropriate to use the 0-H bending frequency (5 x lO 1 ^ Hz) 
for v. The time-scale for transformation observed in the laboratory, 

45 min at 140 K (Hagen, Tielens, and Greenberg 1981), translates then 
into an activation energy of about 5500 K. This is, of course, somewhat 
less than the binding energy of an molecule in an H 2 O ice (-8000 K) 
and actually very similar to the heat of formation of the crystalline 
state from the amorphous state. At a temperature of 90 K, appropriate 
for late-type giants, the transformation time-scale is then about 
2 x 10^ yr. This is much longer than the time-scale for the growth of 
the ice mantle. It is also much longer than the cooling time-scale of 
the grains owing to the outflow (-10^ yr). The ice will condense in the 
amorphous form and will not be able to transform to the crystalline 
state. Thus, although the amorphous state is metastable with respect to 
the crystalline state, the latter is kinetically inaccessible. It is, 
thus, not surprising that the 3 um H 2 O ice band detected in the spectrum 
of the late-type giant 0H0739 has a shape characteristic for amorphous 
H 2 O ice at about 80 K (Hagen, Tielens, and Greenberg 1983b). 

2.5 Polycrystalline Carbon Grains 

A polycrystalline solid consists of a large number of microcrystals 
(10-100 k) separated by some connective material. Within the microcrys- 
tals the molecules are positioned on a regular lattice as in a crystal- 
line solid. The connective material at the grain boundaries is highly 
disordered (Fig. 4). The microcrystalline description of amorphous 
soiids did have some popularity in the past, in particular for glasses, 
but has now generally been deserted for the noncrystalline description 
shown in Fig. 1. The two descriptions differ, of course, in the extent 
of the local order. It should be emphasized that fine-grained polycrys- 
talline solids certainly exist in nature, but experiments have shown 
that in general amorphous solids are not microcrystalline (Zallen 1983; 
Phillips 1981). As will be discussed below, interstellar grains may 
have a polycrystalline structure because of their formation history.. 

Figure 5 shows schematically the structure of a carbon (soot) 
particle. This is essentially a polycrystalline structure. On a small 
scale (diameter < 50 k) the carbon atoms are arranged on a planar honey- 
comb lattice, as is graphite. The carbon atoms are bonded by regular 
covalent sp^ bonds to three neighbors. These bonds are localized. The 
fourth electron of each carbon atom is delocalized among the n bonds 
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Figure 4. Structure of a polycrystalline material. Highly ordered 
microcrystals are separated by highly disordered connective material. 
Note that unlike the amorphous structure shown in Figure 1, there are 
large regions of fairly ordered material present. 


and resonates between several valence-bond structures, giving each C-C 
bond a one- third double-bond character. These small planar "crystals" 
are in fact just the carbon skeleton of polycyclic aromatic hydrocarbons 
(PAHs). Typically, they are less than 50 A in diameter. Several of 
such PAHs are stacked together in platelets bonded by weak van der Waals 
forces, as are the sheets in graphite. Often there is also a consider- 
able content of disorganized, tetrahedrally bonded, sp^ carbon (i.e., 
diamond structure), which distort the planar structure and can link 
different layers within a platelet. Functional groups at the edge of 
these platelets can also cross-link different layers within a plate- 
let. It is these functional groups that separate different platelets 
and prevent orientation of the platelets with respect to each other in a 
soot particle. Although these platelets are topologically similar to 
graphite, they will have quite different optical properties (cf. 
Allamandola, Tielens, and Barker, elsewhere in this volume). These 
platelets are assembled into a much larger, disordered unit, the soot 
particle, bonded by weak van der Waals forces or by bridging aliphatic 
hydrocarbon groups (i.e., tetrahedrally bound carbon skeletons rather 
than aromatic skeletons). As in an amorphous solid, there is thus an 
absence of long-range order, although the structural units themselves 
have highly regular, planar hexagonal structures. 
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Figure 5. Structure of a soot particle. The atoms are arranged on a 
planar honeycomb lattice. Each layer can be viewed as a small graphite 
sheet or the carbon skeleton of a large polycyclic aromatic hydrocarbon 
(PAH) molecule. The structure of benzene, the smallest PAH, is shown in 
detail. Several of these layers form a platelet. These platelets are 
the building blocks of soot particles. 


It should be emphasized that carbon rapidly condensed on a sub- 
strate will form an amorphous solid with aromatic and tetrahedrally 
bound carbon atoms "randomly" interspersed. The differences in bond 
length and angle will produce a disordered structure similar to that 
shown in Figure 1. The soot particles described above are from a struc- 
tural point of view intermediate between this "truly" amorphous carbon 
and graphite. 

These two different views of interstellar amorphous carbon grains 
merely reflect a difference in the condensation history. The polycrys- 
talline particles of Figure 5 will form when many condensation nuclei 
(the PAHs) are formed simultaneously in a (slowly) condensing gas. Most 
of the carbon in the gas will crystallize out into these nuclei. Far- 
ther downstream in the outflowing wind these PAHs can then coagulate, or 
cluster, into the much larger soot particles. In contrast, truly amor- 
phous carbon can be formed when only a few condensation nuclei (i.e., 
PAHs) are formed in a rapidly condensing gas. The condensation of the 
remainder of the carbon forms large grains with an amorphous struc- 
ture. This discussion is, of course, very general, since it is likely 
that the condensing species are molecules or ions rather than carbon 
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atoms. The high abundance of PAHs with about 20 carbon atoms in plane- 
tary nebulae suggests that the first process dominates in that case 
( Allamandola, Tielens, and Barker 1986). The presence of only weak IR 
emission features in the population I WC 9 Wolf-Rayets stars, suggests 
that the second process may actually dominate in these outflows (Cohen 
et al. 1986). In any case, the differences between amorphous carbon and 
polycrystalline carbon are too subtle to be of importance at the present 
state of our knowledge of interstellar carbon grains. However, as the 
observations and theory progress, we anticipate that these differences 
will be revealed. 

2.6 Macromolecular Organic Refractory Grains 

The organic refractory component of the dust in the diffuse interstellar 
medium contains methyl (CHg), methylene (CH 2 ), carbonyl (C=0) and proba- 
bly alcohol (OH) groups. The relative abundance of these groups sug- 
gests that the molecules have a highly branched carbon skeleton or are 
fairly small (Tielens et al. 1986b). The latter possibility is, how- 
ever, ruled out by the requirement that they have to be able to survive 
in the harsh environment of the diffuse interstellar medium. The molec- 
ular structure suggested by the IR data is that of a copolymer (cf. 

Fig. 6). The side groups in this chain can be a hydrogen atom or a CH^, 
C=0, or OH group in random order. Also an occasional backbone C atom 
could be in the carbonyl (C=0) form. Typically, a polymer contains 
about 10^ backbone C atoms but the polymer length is not well known in 
the interstellar case. 
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Figure 6. Monomer repeat unit within a polymer chain. The carbon atoms 
in these units form a linear chain; Y indicates a side group. 
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The structure of such a solid polymer will be very disordered. 
Unlike crystalline solids or monoatomic solids there is no repeatable 
sub-unit (the unit cell). Figure 7 represents a perspective sketch of a 
small section of a polystyrene polymer to give an impression of the 
spatial disorder. In an isolated polymer all the bond lengths and bond 
angles are equal to the crystalline value. Even the third nearest 
neighbors have the correct crystal separation (i.e., staggered bond 
configuration). This set of crystal-like requirements is far stricter 
than for the amorphous structures discussed before (cf. Fig. 1). Yet, a 
polymer has a great deal of configurational freedom, which has been used 
here to open up the polymer as much as possible to show its molecular 
structure clearly. Obviously, a polymer is quite flexible and can 
change directions very often along its chain by rotation about a back- 
bone bond. In that sense, Figure 7 is somewhat misleading since it may 
give the incorrect impression that an amorphous polymer solid consists 
of aligned strands of polymer molecules. In fact, such a solid has a 
random-coil structure, as illustrated in Figure 8. Each polymer adopts 
a random-walk configuration and the structure of the organic glass 
consists, thus, of many intertwined random-coils. Note that the scale 
of Figure 8 is about 1000 times larger than that of Figure 7. The 
bonding between different polymer strands within an amorphous polymer 
solid is due to weak van der Waals forces or hydrogen bonding forces 
between side groups. It should be emphasized that we expect that 
because of their chemical formation process, the polymers in the organic 
refractory component of the interstellar dust are not linear, one- 
dimensional molecules, but rather very branched, three-dimensional 
molecules. That is, they are cross-linked by true chemical bonds 
between side groups (e.g., resins). For that reason they are more 
properly called macromolecules. 


3. OPTICAL PROPERTIES OF AMORPHOUS MATERIALS 


The amorphous nature of materials has considerable influence on their 
physical properties, such as specific heat and thermal conductivity, at 
very low temperatures (Phillips 1981). Essentially, this is due to the 
existence of more than one amorphous state very close in energy. Tun- 
neling between these states can then be of importance at very low tem- 
peratures. Generally, this is not very important at interstellar tem- 
peratures (>10 K). However, the disordered structure of a solid is of 
considerable importance for its optical properties, even at higher 
temperatures. As mentioned in the introduction to §2, it is generally 
thought that a large fraction of the interstellar dust has an amorphous 
structure. This is based on the shape of the 10 um feature and on the 
far-infrared absorption law (-X -1 ). The high near-infrared absorption 
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^ SIDE GROUPS 

Figure 7. Small section of a polystyrene polymer. The carbon atoms are 
tetrahedrally bonded. Alternating carbon atoms in this polymer contain 
a side group (a benzene ring in this case). For clarity, the hydrogen 
atoms have been left out. 



Figure 8. Random coil network for organic glasses. One polymer chain 
has been drawn with a thick line for ease of visualization. 


of circumstellar silicates (Jones and Merrill 1976; Bedijn 1977; Schutte 
and Tielens 1985) may also indicate an amorphous structure (Tielens 
T983). In this section the optical properties of amorphous materials is 
discussed with the emphasis on silicates. The optical properties of 
amorphous carbon grains are discussed by Allamandola, Tielens, and 
Barker elsewhere in this volume. Each of these three observational 




indications of the amorphous nature of interstellar grains will be 
discussed in turn. 

3.1 UV and Visible Absorption 

Absorption in the ultraviolet and visible spectral regions corresponds 
to changes in the electron structure of the material, which is dominated 
by the short-range order. In general, provided the chemical bonding is 
the same in crystalline and amorphous materials, the electronic density 
of states will also be very similar and therefore physical characteris- 
tics such as the energy gap will be about the same. Thus, the fundamen- 
tal optical absorption edge in a crystalline and amorphous material will 
occur at a very similar frequency. For example, the sharp rise in the 
reflection spectra of amorphous and crystalline quartz at about 10 eV 
arises from the onset of the bonding*antibonding transitions of the 
electrons (Tauc 1972). Often, crystalline solids have special 
directions associated with the crystallographic axes, leading to aniso- 
tropy in the optical properties (i.e., dependence on crystal orienta- 
tion). Amorphous materials lack such directionality, that is, they are 
isotropic. 

Crystalline materials can show fine structure in their UV spectrum, 
particularly at low temperatures. Because each unit cell has an identi- 
cal environment, such optical selectivity will arise because energy, as 
well as momentum, has to be conserved in the photon absorption process 
(cf. Kittell 1976). Such spectral structure is an intrinsic property of 
crystalline materials only, because it is a consequence of the long- 
range order (e.g., crystal symmetry). It will be smoothed out in amor- 
phous materials, because of the absence of periodicity (e.g., the wave 
vector is no longer a good quantum number). For example, crystalline 
materials show a rather abrupt onset of the electronic absorption edge, 
which in amorphous materials is much more gradual. Of course, spectral 
substructure in the UV can also result from the presence of more elec- 
tronic transitions within the molecular sub-unit. Such substructure 
will be present in both crystalline and amorphous materials. 

Silicates formed in the outflow from Miras show considerable 
absorption around 1 um (Jones and Merrill 1976; Bedijn 1977; Schutte and 
Helens 1985). This high, near-infrared absorption is somewhat surpris- 
ing since it occurs in the energy gap region; it may indicate an amor- 
phous structure (Tielens 1983a). For an ideal crystalline dielectric 
material, there is an energy gap between the valence band and the con- 
duction band. This energy difference is determined by the bonding- 
antibonding splitting of the electronic states associated with the 
chemical bonding of the nearest neighbors. This energy gap is com- 
pletely devoid of electronic states, and the crystalline material will 
be transparent in this region. An amorphous material will have 
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localized electronic states (e.g. f associated with a few atoms) in the 
energy gap region, and the density of electronic states is, thus, not 
zero but merely low. The transition probability of these states is, 
however, quite small and an amorphous material can be quite transparent 
(cf. glass). As a result, an amorphous material can still be character- 
ized by an optical band gap energy. These localized states are there- 
fore an unlikely source of the enhanced absorption of circumstellar 
silicates at about 1 pm; a more likely source is the presence of impuri- 
ties (e.g., FeO) that can also produce localized electronic states in 
the energy gap region. Unfortunately, there are presently no laboratory 
studies that directly address this astrophysical problem. 

3.2 Mid-Infrared Absorption 

For low frequencies, the oscillations in the electromagnetic field are 
slow enough that molecular vibrations can respond to them, leading to 
absorption. Since the "force constants" between atoms and between 
atomds and electrons are very similar (i.e., both follow from Coulomb 
interaction), it follows from Hooke's law that while electronic transi- 
tions fall in the UV and visible, transitions involving atomic vibra- 
tions will occur in the infrared (i.e., (m e /M) 1/2 - ICT^). Of course, 
the difference in mass also implies that the absorption strength of 
atomic vibrations is smaller by a factor of 10 “^ (consider a harmonic 
oscillator) than that of electronic vibrations (e.g., f IR - 10 _2 -10"3 
whereas f uv - 1-10). Finally, it should be noted that damping in a 
solid is due to collisions with neighboring atoms. The collision rate 
is then approximately equal to the natural frequency of the vibration. 
The damping coefficient is proportional to this collision rate, to the 
anharmonicity of the interaction potential and, of course, to the number 
of phonons present (i.e., temperature). 

In a crystalline solid the vibrational excitations are associated 
with plane waves, and they are therefore characterized by a frequency 
and a wave vector. These lattice vibrations are called phonons. Momen- 
tum conservation is more important for lattice vibrations than for 
electronic excitations, and only a small number of all the possible 
lattice vibrations are optically active. The infrared spectrum of a 
crystalline material consists, therefore, of a few sharp lines. 

The disorder in an amorphous solid localizes the lattice vibrations 
and they can no longer be viewed as plane waves. As a consequence of 
the reduced symmetry, most modes will be infrared active and the infra- 
red absorption spectrum reflects to a greater extent than it does in 
crystals the density of vibrational states. It should be emphasized 
that the vibrational density of states in an amorphous material is not 
much different from that of the corresponding crystalline material. 
Basically, it reflects the interaction between nearest neighbors, which 
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is very similar in both types of materials. Of course, the increased 
dispersion in bond lengths and bond angles will produce a larger range 
in bond strength and thus in frequency (recall Hooke's law). Because of 
this and because most modes are IR active, the fundamental absorption 
features in an amorphous material will be broader than in a crystalline 
material. 

The integrated strength of an absorption band in an amorphous 
material is generally somewhat less than that of the crystalline mate- 
rial. That the oscillator strength is very similar is of course not 
surprising, since the force constants of the oscillators are very simi- 
lar. The larger disorder in an amorphous material implies, however, a 
larger damping coefficient (i.e., a larger effective coupling with 
phonons). As a consequence the peak strength of the absorption in an 
amorphous material is smaller and the width larger (cf. a damped oscil- 
lator). The absorption strength is related to the change in dipole 
moment. The high order in 'a crystalline material makes all dipoles 
change in unison, resulting in a fairly large change in dipole moment. 

In contrast, the disorder in an amorphous material will make the effec- 
tive change in dipole moment less and consequently lower the integrated 
absorption strength (i.e., the oscillator strength). For silicates this 
difference in integrated absorption strength of amorphous and crystal- 
line solids seems to be relatively large (a factor of 2.5; see §1.1). 
This is much larger than in, for example, amorphous ice (25 Hagen 
et al. 1981). It should be noted, however, that the variation in peak 
absorption strength and width are much larger than the variations in the 
integrated strength. In general, it is therefore better to use the 
integrated strength in interpreting the strength of absorption features 
(cf. §1.1). 

Finally, it is interesting to consider the effects of temperature 
on the fundamental vibrations of a solid. Typically, absorption fea- 
tures sharpen and shift to slightly higher frequencies with decreasing 
temperature, but the integrated absorption strength stays approximately 
constant (Plendl 1970). These variations are a result of the thermal 
shrinking of the lattice. The decreased nearest neighbor distances 
imply larger force constants and thus higher frequencies (Hooke's 
law). The decrease in width is due to the decreased importance of 
anharmonicity and to a decrease in the number of phonons present, imply- 
ing a smaller damping coefficient (cf. a harmonic oscillator). Differ- 
ent modes (e.g., stretching versus bending vibrations) in a material may 
behave differently upon cooldown, since they may have different anharmo- 
nicity constants and they may interact with different phonons. This may 
account for the experimentally measured behavior of amorphous magnesium 
silicates (Day 1976). 

At this point it is tempting to speculate that the decreased width 
of the 10 urn feature in absorption in the interstellar medium 
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(presumably at about 10-20 K) as compared to the emission feature in 
Orion (a250 K; Gillett et al. 1975b; Roche and Aitken 1985) is due to 
the effects of thermal shrinking. However, no experimental data are 
available on the width of the 10 ym silicate band at low (10 K) tempera- 
tures. The variation seems, however, consistent with that measured for 
other materials (Plendl 1970). 

3.3 Far-Infrared Absorption 

The far-infrared (X > 30 ym) absorption of amorphous materials differs 
considerably from that of crystalline materials. In crystalline mate- 
rials most of the far-infrared absorption is due to multiphonon pro- 
cesses (for example, one photon and one phonon in, and one phonon of 
higher energy out). Fundamental vibrations (i.e., one photon in, one 
phonon out) are optically inactive because of wave-vector conservation 
(cf. §3-2). Multiphonon processes are highly temperature sensitive, and 
at 10 K the absorptivity will be much reduced, essentially because there 
are few phonons around (Hadni 1970). Besides multiphonon processes, a 
material can also absorb in the damping wing of an infrared active 
fundamental vibration. Far from the natural frequency, the absorption 
coefficient of a damped harmonic oscillator will show a quadratic depen- 
dence with frequency, and such a dependence is expected for crystalline 
dielectric materials at low temperatures (Wooten 1972). In a metallic 
or semimetallic material (e.g., graphite), far-infrared absorption is 
due to interaction of the radiation with the free electrons (Wooten 
1972). Again, the absorption coefficient will depend quadratically on 
the frequency, since the free electrons in a material are essentially 
"oscillators” without a restoring force. The damping is now due to 
scattering of the electrons by phonons (i.e., lattice vibrations). 

In amorphous materials, the far infrared absorption is quite 
strong, often featureless, and not temperature sensitive (Hadni 1970). 
These differences with respect to crystalline materials are due to the 
loss of long-range order. The resulting breakdown of the selection 
rules implies that all modes are infrared active, and the far-infrared 
absorption spectrum reflects the density of vibrational states. Gener- 
ally, the density of vibrational states is well described by the Debye 
law at low energies (~w 2 ; Kittell 1976). The strength of the far- 
infrared absorption depends on the disorder in the charge distribution 
in the lattice. But, since the oscillator strength is independent of 
the particular vibrational mode involved in the phonon regions, the 
frequency dependence of the far-infrared absorption reflects directly 
the phonon spectrum (Tielens 1986a). Thus, an amorphous material will 
show the same wavelength dependence of the far-infrared absorption 
(i.e., x" 2 ) as a crystalline material, albeit for completely different 
physical reasons. The far-infrared absorption of amorphous materials 
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is, of course, much stronger and not temperature sensitive. The far- 
infrared spectrum of fused quartz shows these characteristics of amor- 
phous materials (Hasegawa and Koike 1984). 

An exception to the quadratic law for far-infrared absorption 
should be made for amorphous, layered materials, such as amorphous 
carbon, and layer-lattice silicates. The structure of such materials 
limits the phonons to two dimensions and, consequently, their phonon 
spectrum is proportional to the frequency rather than to the frequency 
squared (Kittell 1976). The temperature dependence of the specific heat 
of graphite, the classic example of a layered structure, reveals the 
presence of such a X " 1 phonon spectrum (Touloukian and Buyco 1970). 
Because the far-infrared absorption of this material is dominated by 
free electrons, the far-infrared spectrum of this crystalline, layered 
material shows a X " 2 dependence (see above). But, amorphous carbon, 
which has a very similar structure (cf. §2.5) does show the expected 
x" 1 far-infrared absorption law (Koike et al. 1980). Layer-lattice 
silicates, which also have a layered structure, show a X“ n far- 
infrared absorption law with n in the range 1.25 to 1.5 (Day 1976). 
Apparently, there is some degree of cross-linking between the layers in 
such materials, leading to a somewhat steeper law than expected for a 
layered structure (Tielens 1986a). The last case of a X -1 far- 
infrared absorption law involves small amorphous grains. In small 
particles (say <30 A) the surface modes (two-dimensional) rather than 
the bulk modes (three-dimensional) will dominate the far-infrared 
absorption (Seki and Yamamoto 1980). Again, for amorphous materials 
where all modes are active, the two-dimensional structure will lead to 
a X"* 1 law. 

Clearly, the far-infrared is a very useful tool for studying the 
structure of interstellar grains. Far-infrared observations of celes- 
tial objects generally show a far-infrared absorption law that is less 
steep than expected fo.r crystalline materials (Campbell et al. 1976; 
Gatley et al. 1977; Harvey et al. 1978). This implies that most of the 
interstellar grains have an amorphous layered structure. Since amor- 
phous carbon seems to be a relatively minor component of the inter- 
stellar dust (Table I), this suggests that interstellar silicates, which 
may make up half of the interstellar dust volume, have a layer-lattice 
amorphous structure. In this respect, it is worth noting that the 
2.9 urn band observed in the spectrum of the galactic center can be well 
fitted by a spectrum of hydrated magnesium silicates, a layer-lattice 
silicate (Tielens et al. 1986b). Obviously, a search for this band in 
other sources with strong silicate bands is important for the corrobora- 
tion of this conclusion. 
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4. PHYSICAL PRINCIPLES OF GRAIN-SURFACE CHEMISTRY 


Figure 9 shows a schematic view of grain-surface chemistry. Atoms and 
molecules are accreted from the gas phase. They may diffuse and react 
on the surface, and the reaction products may remain on the surface, 
forming a grain mantle, or may actually be ejected from the surface by 
some means. In this section we will briefly discuss some of the physics 
of these processes. For a detailed discussion, the reader is referred 
to the tutorial by Watson (1976). This discussion will be geared toward 
Hj formation and the formation of icy grain mantles in molecular clouds. 


ACCRETION 
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REACTION 


OOOOOOOOOOOOOOOOOOOOOO EJECTION 

Figure 9- Grain-surface chemistry. Gas-phase species accrete, diffuse, 
and react on a grain surface. Possibly, the reaction product is ejected 
into the gas phase. 


4.1 Surface Structure 

The structure of the surface is of importance for many aspects of grain- 
surface chemistry. It should be realized that the structure of the 
surface will deviate from that of the bulk, in particular for a crystal- 
line solid. In a sense, the surface is adjacent to a row of lattice 
vacancies, and as a result the equilibrium positions of the molecules at 
the surface will have been relaxed from that in the bulk. Moreover, the 
surface will facilitate the formation of lattice defects, which will 
further disturb the molecular bonds. Furthermore, the surface composi- 
tion may differ from that of the bulk. Finally, a larger concentration 
of impurities may be present on the surface than in the bulk. Some of 
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these points are illustrated in Figures 10 and 11, which show some of 
the structures possible at silicate and carbon surfaces. 
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Figure 10. Some of the structures possible at silicate surfaces. Note 
the presence of silanol (Si-OH) and siloxane (Si-O-Si) groups. The 
former bonds 1^0 by hydrogen bonding. 
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Figure 1 1 . Some of the structures possible at the edge of the basal 
plane in graphite or amorphous carbon. 
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Surface groups on silicates are extensively discussed in the mono- 
graph by Hair (1967) and by Snoeyink and Weber (1972). In an astrophys- 
ical context they have been briefly discussed by Duley and McCullough 
(1977). Generally, a silicate surface contains two types of surface 
groups, silanol groups (Si-OH) and siloxane groups (Si-O-Si). Silicates 
precipitated from aqueous solutions at low temperature (<700 K) mainly 
contain hydroxyl groups on the surface, a result of the tendency of the 
silicon atom to complete tetrahedral coordination. Water will be 
adsorbed on these surface silanol groups by van der Waals and hydrogen 
bonding forces. Upon heating to 700 K, these physically adsorbed HgO 
molecules will be lost (dehydration of the silicate surface; Fig. 10). 
Adjacent hydroxyl groups will also form H 2 O molecules that are ejected 
from the surface, leaving a strained siloxane group behind. Because of 
the strain in these siloxane groups, rehydrogenation can occur, reform- 
ing the silanol groups. Warm-up above 700 K (or formation above this 
temperature) enables the strained siloxane groups to anneal themselves 
and, upon addition of water, no rehydrogenation will occur (e.g., the 
dehydrogenation is irreversible). Above 1000 K, only isolated silanol 
groups (so-called free OH or terminal OH groups) and siloxane groups 
remain (Fig. 10). Such a surface will bind H 2 O much more weakly than 
the untreated silicates and is actually quite hydrophobic. The presence 
of other atoms in the silicate surface structure, such as Al, will 
increase the acidic nature of a silicate surface because of the induc- 
tive effect on adjacent silanol groups. Such atoms may also bind 
hydroxyl groups, which will behave slightly differently on heat 
treatment. 

Physical adsorption of a polar molecule, such as HjO, on a silicate 
surface will preferentially occur on the OH groups, and the heat of 
adsorption will increase with increasing polarity of the adsorbed mole- 
cule, a result of hydrogen bonding. The silanol groups on a silicate 
surface react as a weak acid and can, thus, react with bases (at room 
temperature). Reactions with alcohols (at high temperatures) will 
produce estersils (e.g., Si-0-C groups). At low temperatures, such 
reactions are inhibited by activation energies. Annealed siloxane 
groups are very weak bases and show, therefore, only a weak tendency to 
form hydrogen bonds. They are also quite unreactive. Because inter- 
stellar silicates condense out at relatively high temperatures 
(-1000 K), they are presumably initially quite hydrophobic. After 
processing in a shock they will presumably have many dangling surface 
bonds, which will be quickly saturated by H atoms. The resulting sur- 
face may bind H 2 O more strongly. In any case, the low temperature of 
interstellar dust grains will inhibit reactions with surface molecules, 
and the only property of interest in interstellar grain surface chemis- 
try is the ability of silicate surfaces to bond H 2 0 and other polar 
molecules relatively strongly by hydrogen bonding. 
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Commercial carbon blacks usually have several different surface 
groups besides aromatic hydrogen (e.g., methyl, phenol, and carboxyl 
groups; Fig. 11). These groups are substituted for the aromatic hydro- 
gen at the edge of the aromatic ring system (Puri 197.6). They result 
from chemical reactions with the graphitic structure (at elevated cem- 
peratures). Graphite or amorphous carbon freshly injected into the 
interstellar medium will contain mainly aromatic hydrogen at its edges 
because of the overabundance of hydrogen in their birth sites. Sputter- 
ing in fast shocks will change this surface structure and may lead to 
dangling bonds. These will either be resaturated by atomic hydrogen or 
the carbon atoms in the edge may revert to a divalent s 2 p 2 state. The 
latter structure limits, of course, the resonance of the ir electrons 
and makes that site more reactive, but leaves no dangling bonds. In 
fact, vibrational spectroscopy of fluorescing interstellar PAHs, the 
building blocks of the carbon grains, generally shows an absence of 
surface groups other than aromatic hydrogen ( Allamandola, Tielens, and 
Barker, elsewhere in this volume), and this is likely to be true for 
carbon grains as well. 

4.2 Adsorption Energies 

A species can be either physically or chemically bonded to the sur- 
face. In physical adsorption the interaction is due to van der Waals- 
London or electric field-dipole forces, and the interaction energy is 
fairly low (<0.5 eV). There is no charge exchange between the adsorbed 
species and the substrate in physical adsorption. The attractive force 
is due solely to the interaction of the (instantaneous) dipole of the 
adsorbed species and the neighboring surface molecules (Kittell 1976). 

In chemisorption, the interaction leads to actual chemical bond forma- 
tion and the interaction energy is much larger (-0.5-5 eV). This inter- 
action can entail a complete charge exchange between the surface and the 
adatom and the formation of an ionic bond. Usually, however, the elec- 
tronic wave functions of the adatom and the surface mix into a new wave 
function. Essentially, the electrons responsible for the bonding move 
in orbitals between the surface and the adatom, and a covalent bond is 
formed . 

The theoretical description of chemisorption is very complex and 
not completely understood (cf. Einstein et al. 1980). For illustration 
purposes consider the physical processes occurring when an adatom 
approaches a metal surface. A crude impression of the tendency to form 
a chemical bond can be obtained by comparing the adsorbate's ionization 
potential and electron affinity with the work function of the surface. 

If the ionization potential of the adatom is less than the work function 
of the surface, then an electron can be transferred to the surface and 
the adatom becomes positively charged. Charge transfer can also take 
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place when the work function is less than the electron affinity of the 
adatom. Now the electron is transferred to the adatom. In both cases, 
chemisorption is energetically feasible. However, if the work function 
of the surface is larger than the electron affinity but less than the 
ionization potential of the adatom, then the adatom will be physically 
adsorbed. This discussion is of course very crude because the interac- 
tion between an approaching atom and a surface will lead to broadening 
of the atomic states, and the ionization and affinity levels will become 
bands whose widths depend on the degree of interaction (cf. Einstein 
et al. 1980). These bands may partially overlap with the (Fermi) level 
of the metal. When the atom is very close to the surface, strong inter- 
action may actually lead to a complete modification of the electronic 
structure of the atom and the surface. 

A simple potential energy diagram for chemisorption is shown in 
Figure 12. Note that the physical adsorption forces have a longer range 
than the chemical forces and, thus, an accreting species will first 
encounter the physical adsorption potential. For atoms approaching a 
surface there is no activation barrier associated with chemisorption. 

For molecules, however, there is probably an activation barrier because 
a molecular bond has to be weakened or broken. In the interstellar 
medium, chemisorption of molecules is, therefore, generally unimportant 
owing to the low temperatures. 




Figure 12. A potential energy diagram for the adsorption of an atom (a) 
or a molecule (b) on a surface. Note that there are two adsorption 
wells, the physical adsorption well with binding energy, Up and the 
chemisorption well with binding energy, U c . For molecules, these wells 
may be separated by an activation energy, E. 
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to 2 x lO^/cra^, corresponding to 


Hydrogen atoms on silicate surfaces have two types of binding 
sites. First, a chemisorbed site with a binding energy of about 2 eV 
(-2.3 * 10^ K). Second, a physical adsorbed site with a binding energy 
of only about 0.09 eV (-1000 K; Gelb and Kim 1971; Wood and Wise 1962; 
Hickmott I960; de Boer and van Steenis 1952; Langmuir 1912, 1915). The 
second sites should be viewed to be on top of the first one (cf. 

Fig. 12) and will only be occupied when the first one is already 
filled. Mote that reaction between H atoms in the two layers is inhib- 
ited by an activation barrier (-1000 K) associated with breaking the 
chemisorbed bond with the surface (see §5). The number of sites on a 
silicate surface is about 1 x 10^ 
about 7 A 2 per site. 

Although it is well known that atomic hydrogen can be chemisorbed 
on the basal plane of graphite as well as on the edges (King and Wise 
1963; Wood and Wise 1969; Rye 1977; Gould 1975; Balooch and Olander 
1975; Robell et al. 1964), the parameters of these states are less well 
known, a result of the complicated chemistry of H atoms with the edge 
structure of such surfaces at elevated temperatures (300 K). At high 
surface coverage, H 2 formation has an activation barrier of about 
1200 K, suggesting an H chemisorption energy on the basal plane of 
perhaps 2 eV. From diffusion experiments, McCarroll and McKee (1971) 
estimated a somewhat lower binding energy of about 0.9 eV on the basal 
plane. For H 2 formation on interstellar grains, the activation energy 
is, however, more important than the binding energy (see §5), so that 
this discrepancy should not worry us too much. The "aromatic" H atoms 
at the edge of the basal plane are chemically bonded by about 5 eV. The 
physical adsorption energy of H atoms on the basal plane of graphite is 
measured to be about 800 K (Lee 1975). The number of chemisorption 
sites is again of the order of 10 3 cm . Other atoms on silicate and 
graphite surfaces are also expected to be chemisorbed on these sites. 

At low temperatures, molecules will be physically adsorbed on these 
surfaces . 

Atoms and molecules will be physically adsorbed on icy grain. 

Their binding energies can be determined from molecular beam studies and 
from diffusion experiments in low-temperature matrices. Table II summa- 
rizes estimates for the binding energy of relevant atoms and molecules 
on an icy surface. 

Using these binding energies, the residence time of a species on a 
surface can be calculated 


t = t q exp(E b /kT) , (4) 

where t q is the oscillation time of the adsorbed species perpendicular 
to the surface (Frenkel 1924). Essentially, the adsorbed molecule has 
replaced at least one degree of translationary freedom by a vibrational 
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TABLE 2 


The Interaction of Atoms and Molecules with an H 2 0 Surface 


Species 

Eb(K) 1 

t ev (sec) 2 
10 30 

5 

x m (sec) 3 

10 

30 

Notes 

H 

350 

1.6*10 3 

10" 7 

10~ 12 

io- 12 

10” 12 

4 

H 2 

450 

3*10 7 

3*10" 6 

4x10' 12 

4x10“ 12 

4x10~ 12 

4 

C 

800 


4x10 -1 

10 8 

10" 2 

2x10“ 9 

5,6 

N 

800 


4x10“^ 

10 8 

10“ 2 

O' 

j 

o 

X 

CM 

5,6 

0 

800 


4x10 _1 

10 8 

10~ 2 

2x10* 9 

5,6 

S 

1100 


8x10 13 


10 2 

4x 10~ 8 

6,7 

CO 

1900 




6xio 12 

2x10* 4 

5,8 

N 2 

1700 




10 10 

2x10‘ 5 

5,8 

°2 

1600 




7x10 8 

9x10 -6 

5,8 

ch 4 

2600 





2X10” 1 

5,8 

h 2 o 

4000 





2x10 5 

9,6 


Notes: 1) Estimated binding energy of a species on an H 2 0 surface. 

2) Evaporation time-scale from an H 2 0 surface for grain 
temperatures of 10 and 30 K. No value is given when it exceeds the 
lifetime of a molecular cloud. 

3) Migration time-scale on an H 2 0 surface for grain tempera- 
tures of 5, 10 and 30 K. No value is given when it exceeds the lifetime 
of a molecular cloud. 

4) Binding energy determined from experimentally measured 
attractive potential corrected for zero point energy and H 2 coverage 
(see text) . 

5) Binding energy determined from measured maximum attractive 
potential for H 2 , corrected for difference in polarizability (cf. Watson 
1976). 

6) Diffusion barrier determined from low temperature matrix 
isolation studies (Tielens and Hagen, unpublished). 

7) Binding energy determined from diffusion experiments in 
low temperature matrices in analogy to oxygen. 

8) The diffusion barrier is assumed to be 30% of the binding 

energy . 

9) Enhanced binding energy due to hydrogen bonding with the 
H 2 0 surface. 
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degree of freedom of the bond with the surface. It may also have lost 
some degrees of rotational freedom. For a monoatomic gas which is 
completely mobile on the surface, we can find for t q from a statisti- 
cal thermodynamical analysis 


, 0 = (h/KT)f z 


(5) 


where f z is the vibrational partition function of the adsorbed mole- 
cule (de Boer 1968). At high temperatures the vibration perpendicular 
to the surface is excited and f z is equal to kT/hv z , where v z ^s 
the vibrational frequency perpendicular to the surface and r = v“ . 

We will assume that the frequency of vibration perpendicular to the 
surface is equal to that parallel to the surface, which for a symmetric 
harmonic potential is given by (Landau and Lifshitz I960): 


v 2 = [(2N s E b )/ir 2 m] 1/2 


( 6 ) 


where N g is the number of sites per unit surface area and m the mass 

of the atom. With a physical adsorption energy of 1000 K and 

N * (7 A 2 )” 1 , v_ is about 5 * 10 l2 sec" 1 for an H atom and about 
s z 

10 12 sec" 1 for an 0 atom. For a chemisorbed site with a binding energy 


of 2 eV, v z is about 2 x lO 1 ^ and 5’* 10 l<: see - ' , respectively. At low 
temperatures, the vibration perpendicular to the surface is not excited 
and f z = 1, leading to t q = h/kT • 5 * 10“ 2 sec, where a temperature 
of 10 K has been assumed. When the adsorbed species is not mobile on 
the surface or has lost some rotational degrees of freedom, then equa- 
tion (5) has to be multiplied by the ratio of the partition functions 
for hindered and free translationary motion on the surface or the rota- 
tional partition functions on the surface and in the gas phase, respec- 
tively. In both these cases t q will be less than h/kT. In the limit 
that an adsorbed atom has lost all its degrees of translationary free- 
dom, a different expression is derived: 


12 


.-1 


t = ( h/kT ) f „f „f _ ( N h /2itmkT ) 

O a V 2 3 


(7) 


where the f x , fy are the partition functions along the surface 
(de Boer 1968). The last factor in this expression takes into account 
the number of possible ways of distributing the adsorbed species over 
the adsorption sites. Consider again the low-temperature limit 
(T = 10 K), assume a surface area per site of 7 A 2 , and a mass of 
— This yields t * 4 xJO -1 sec. Thus, for the interstellar 


20 amu . 
case, 


lies between 


10 


-13 


and 4 x 10 


-12 


sec for an H atom physi 

cally adsorbed and chemisorbed, respectively. Using t q 
find the residence time-scales on an icy surface given in Table II. 


- 12 

=10 sec, we 
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4.3 Sticking Probabilities 

A gas-phase species approaching a grain experiences an attractive force 
owing to the molecules in the grain. Whether a gas-phase species will 
stick or not depends on two factors: (1) the magnitude of its transla- 
tional energy compared to the adsorption energy and (2) the efficiency 
with which the excess kinetic energy can be transferred to the surface 
atoms with which it is colliding. Some insight into the physics can be 
gained by examining the classical case. Imagine an atom with mass m 
colliding with a single surface atom of mass M which is attached to 
the surface with a spring of a frequency . (the Debye frequency). 

This is the so-called "soft-cube" model. We can define a characteristic 
inverse time-scale for the collision, w Q , 

u) Q = v/2b = [(E th + E b )/2mb 2 ] 1/2 (8) 

where v is the velocity during the approach, b is a characteristic 
length scale for the repulsive interaction (»0.3-0.4 A), and E fch and 
E b are the kinetic and binding energies, respectively. For an H atom 
at low temperatures, <u 0 is typically 3 * 10 ^ and 6 x 10 ^ sec" 1 for 
chemisorption and physical adsorption, respectively. When the collision 
time is short compared to the spring period (i.e., a» Q >> then the 
energy transferred in the collision is 

AE = 4(m/M)(E th + E b ) (9) 

(Hollenbach and Salpeter 1970). This transferred energy should be 
compared to the initial kinetic energy of the impacting atom, E th , to 
determine whether sticking will occur. If the collision time-scale and 
the spring period are equal (i.e., o) Q = w^) then the energy transferred 
in the collision is one fourth that given by equation (9). If u> 0 is 
actually much less than the spring frequency, then other surface atoms 
can couple in the collision. The energy transferred is then decreased 
over that given by equation (9) by a factor of about 1/2 (u> q /<d d )^ 
(Hollenbach and Salpeter 1970; Watson 1976). Consider, for example, an 
H 2 molecule with a kinetic energy corresponding to 300 K colliding with 
an ice surface (E b » 800 K) . The maximum energy transfer is 50% of the 
total energy (cf. eq. [9]) or somewhat less than 600 K. Actually u> D 
is about equal to w Q in this case, and the energy transferred is only 
about 140 K, leaving about 160 K of kinetic energy, and sticking will 
not occur in the initial collision. Because the collision is generally 
not head-on, the approaching atom will, however, ricochet off and col- 
lide with other surface atoms before being finally ejected from the 
surface; therefore, more energy can be transferred (Hollenbach and 
Salpeter 1970). Taking this into account, the sticking coefficient of 
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an H 2 molecule on an H 2 0 surface at 300 K is calculated to be about 0.3 
(Hollenbach and Salpeter 1970), close to the experimental value (Govers 
et al. 1980). Of course when the surface is covered with a layer of K 2 
molecules, as is likely inside molecular clouds, then the sticking 
coefficient will increase to about unity (i.e., AE * E fch E^), in 
agreement with experiments (Govers et al. 1980). 

In general, if the kinetic energy of the accreting atom is much 
less than the adsorption energy, then only a relatively small amount of 
energy needs to be transferred to the surface, and the sticking coeffi- 
cient will be relatively large. Similarly, if the accreting atom is 
much heavier than the surface atoms, then the amount of kinetic energy 
transferred to the surface will be relatively large and the sticking 
coefficient will also be relatively large. Finally, for a stiffer 
lattice (e.g., u» D >> ai Q ) the sticking coefficient will be lower than for 
a softer lattice. It should be noted that this classical description 
assumes that any amount of energy can be transferred in the collision 
and, thus, when the initial kinetic energy goes to zero, the sticking 
coefficient will go to unity. This is, however, not correct, and for 
low initial energies quantum effects have to be taken into account (cf. 
Leitch-Devlin and Williams 1985). The adatom has to wind up in a bound 
level in the interaction potential well, and, thus, a minimum amount of 
energy has to be transferred, which depends on the spacing of the upper 
levels in the potential well. Because of anharmonicity this will be of 
little importance for chemisorption, but it can be of importance for 
physical adsorption. It should also be kept in mind that the phonon 
spectrum of a solid decreases with energy (e.g., the phonon density is 
proportional to (w/w_) , so that coupling and thus energy transfer to 
the lattice will be less at low impact energies. Again, this is mainly 
important for physical adsorption for which the interaction potential is 
quite small. 

The sticking coefficients can be found from the theoretical and 
experimental results. Anticipating our discussion later on, the stick- 
ing probabilities of importance for interstellar grain chemistry are the 
collision of a species with a clean silicate or graphite surface on 
which it can be chemisorbed (e.g., in the diffuse interstellar medium) 
and the collision of a species with an icy grain mantle on which it will 
be physically adsorbed (e.g., inside dense molecular clouds). Theoreti- 
cal calculations on the sticking probability of an atom on a grain 
surface have been evaluated for several likely gas-phase-species and 
grain combinations (Hollenbach and Salpeter 1970; Watson and Salpeter 
1972; Burke and Hollenbach 1984; Leitch-Devlin and Williams 1985). 

These calculations show that an H atom colliding with a silicate or 
graphite surface on which it can be chemisorbed will have a sticking 
coefficient of unity at a gas temperature of 100 K. An H atom colliding 
with a clean H 2 0 ice surface at 10 K has a relatively low sticking 
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coefficient (*0.3). However, inside a molecular cloud the icy surface 
will be partially covered with H 2 molecules, and this will raise the 
sticking coefficient to close to unity. For heavier atoms, the larger 
binding energy and larger mass will insure a large sticking probability 
( si ) . 

An exception to these high sticking probabilities may have to be 
made for ions accreting on charged grains (Watson and Salpeter 1972). 
Grains inside molecular clouds are negatively charged most of the time 
because the electrons have a greater velocity than the ions (Spitzer 
1978). Charge recombination at or near the surface complicates the 
evaluation of the sticking coefficient considerably (Watson and Salpeter 
1972). The Coulomb interaction energy is much larger than the adsorp- 
tion energy of neutrals on icy grains considered above (i.e., comparable 
to a chemisorption energy). Thus, an accreting ion acquires a much 
larger approach velocity than a neutral. The Coulomb interaction may 
lead to the tunneling of an electron through the work function of the 
grain material, followed by ion-electron recombination, possibly before 
the ion has reached the surface. The resulting neutral experiences the 
much weaker, neutral-grain adsorption potential, but its kinetic energy 
is still that of the original ion. This is, therefore, a case in which 
the kinetic energy of the gas-phase species is much larger than the 
adsorption energy and, consequently, the sticking probability of ions 
will be very low (Watson and Salpeter 1972). This hinges, however, on 
the mobility of electrons in grains, and if electrons are actually quite 
localized in a grain then the ion-sticking probability will be close to 
unity. Of course, the Coulomb interaction will also enhance the colli- 
sion cross section by about an order of magnitude (i.e., focusing). 

A parameter of importance for molecule formation on grain surfaces 
is the time-scale at which a grain accretes species from the gas 
phase. This accretion time-scale is given by 


T ac = ( V vS) 


-1 


( 10 ) 


10- 4 n o ) 


where n^^ is the density of the species in the gas phase (5 
o is the geometrical cross section of the grain (10 - ^ cm 2 ), v the 
average velocity (1 km/sec), and S the sticking probability (*1). 

thus, inside dense molecular clouds a grain will see a very small 
flux of incoming atoms or molecules (e.g., a few per day at a density of 
10 4 cm“3). Of course, with this flux all the heavier species in the gas 
phase will have accreted in 10^ yr, which is a much shorter period than 
the cloud's lifetime. 
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4.4 Surface Mobility 


The binding energy of an adsorbed species varies considerably across the 
surface owing to the discrete nature of the lattice. For physical 
adsorption, the barrier against diffusion to an adjacent site, E^, is 
typically about 30% of the binding energy, although sometimes much lower 
values have been calculated (Jaycock and Parfitt 1986). For chemisorp- 
tion these variations can be comparable to the binding energy, although 
often it is again about half the binding energy. For a heavy species 
the migration time-scale, t m (i.e., the time to go to a neighboring 
site), is given by the Boltzman expression: 

t = v” 1 exp(E d /kT) . (11) 

Here v z is the frequency of the vibrational motion along the surface 
(cf. eq. [6]). For atomic hydrogen, deuterium, and helium and molecular 
hydrogen, migration is due to quantum mechanical tunneling. The time- 
scale for this process is given by 


t = 4h/AE , 
tn 


( 12 ) 


where AE is the width of the lowest energy band (cf. Hollenbach and 
Salpeter 1970; Watson 1976). For atomic H on an 1^0 surface this has 
been evaluated to be 30 K, yielding a surface diffusion time-scale of 
10 -1 ^ sec. Values for the migration time-scale on an. icy surface are 
given in Table II. It should be stressed that the actual time involved 
in moving from one site to a neighboring site is very small 
(» 10“ 12 sec). The limiting step in the migration process is the 
acquisition of the energy required to overcome the diffusion barrier. 

The surface diffusion process is essentially a random-walk one. 

The rms value of the displacement vector R, after a time At is then 


<R 2 > 1/2 


. * 1 / 2 - 
N l 


(13) 


where ft, is the length of one site and N the number of hops during 
this time (N = At/x ). Thus, the average time to scan a grain surface, 
t s , is t = A T m /4n, where A 3 is the total number of sites on the 
surface (-10^ for a 1000 A grain). 

4.5 Surface Reactions 

Here we will only consider the formation of icy grain mantles inside 
dense clouds. The reactions associated with ^ formation are discussed 
in §5. The reaction network can be limited considerably by a comparison 
of the migration time-scale and the thermal evaporation time-scale 



(Table II) with the accretion time-scale (eq. [10]; ~10“* sec). Species 
with a migration time-scale much less than the accretion time-scale can 
diffuse on the surface until they find a co-reactant (e.g., H, H 2 , C, N, 
0, S at 10 K). Otherwise they are trapped in their site and can only 
react with migrating co-reactants that visit their sites (e.g., CO, N 2 , 
0 2 , Clfy, NH^, and H 2 0 at 10 K). Thus, at 10 K the grain mantle reaction 
network contains reactions of the migrating atoms H, C, N, and 0 among 
themselves and with the non-diffusing atoms or radicals (Tielens and 
Hagen 1982). 

Comparison of the evaporation time-3cale with the accretion time- 
scale shows that if atomic H does not find a reaction partner on the 
surface then it will evaporate before the next radical accretes. Direct 
recombination of two H atoms to form H 2 does not occur on icy grain 
surfaces in molecular clouds (cf. §5). It should be stressed that 
because of their low mass, H atoms can tunnel through activation bar- 
riers and react with molecules containing an even number of electrons 
(e.g., non-radicals), such as CO, 0 2 , H 2 C0, H 2 0 2 , H 2 S, and 0^ (Tielens 
and Hagen 1982). The time-scale for tunneling through a barrier with a 
height, E_ , and a width, a, is given by 

a 

T fc = t q exp[(2a/h)(2mE a ) 1/2 ] (14) 

where m is the mass of the tunneling species. Assuming a barrier 
height of 2000 K and a width of 1 A, this time-scale is 2 x 10" 4 sec. 
This is much longer than the migration time-3cale, and the probability 
for reaction is therefore very small (t /T fc * 5 * 10” y ). Thus, gener- 
ally an H atom will migrate to an adjacent site rather than react. 
However, over an evaporation time-scale, an H atom will visit each site 
many times (t /t * 10 U at 10 K for a 1000 A grain). The probability 
of reaction p r is therefore given by 

Pr 1 (,5) 

For the numerical example considered above, p p is about 5. Thus, 

H atoms can react with molecules even if there is an activation barrier 
as large as 2000 K. For a high concentration of co-reactants, p r will 
be larger. For example, for a concentration of \% on a grain of radius 
10^ A, the limiting barrier an H atom can tunnel through is about 
3500 K. Similarly, reactions of with an activation barrier less than 
about 4500 K can be important on grain surfaces (Tielens and Hagen 
1982). Heavier atoms or radicals will have to thermally hop over an 
activation barrier and reactions with activation barriers larger than 
the barrier against diffusion are effectively inhibited unless the 
abundance of co-reactants is very large. Assuming the surface is cov- 
ered with such co-reactants and considering the accretion of a new 


radical as the limiting step for reactions of 0 atoms (rather than 
evaporation), reactions of 0 atoms with an activation energy of about 
450 K or less have to be taken into account at 10 K. Reactions of 
0 atoms with negligible activation barriers include O 2 and possibly CO 
forming 0g and CO 2 , although the latter activation barrier is still 
somewhat controversial (Tielens et al. 1986a). 

4.6 Dissipation of Reaction Heat 

The product molecules of grain-surface reactions are generally formed in 
electronically and vibrationally excited states. The subsequent decay, 
ultimately to the ground state, stabilizes the newly formed molecule. 
Some of the reaction heat might be transferred into center-of-mass 
translational motion, possibly leading to ejection of the molecule. The 
excess translational energy of the ejected molecule will heat the gas. 
The ejected molecule may also contain excess vibrational or rotational 
energy, which will decay radiatively. This process might explain the 
relatively high excitation of J > 4 levels of H 2 observed in diffuse 
interstellar clouds, although other mechanisms are also possible 
(Spitzer and Cochran 1973; Dalgarno 1975). 

In particular, the ejection of "hot" H 2 from grain surfaces has 
incited several theoretical studies. Classical mechanical trajectory 
calculations have been performed, mainly for H 2 formed by recombination 
of two physically adsorbed H atoms on a grain surface (Hollenbach and 
Salpeter 1970; Hunter and Watson 1978). In these calculations enough 
energy was transferred to the center-of-mass translational motion in the 
collision between the newly formed H 2 molecule and the surface potential 
to cause ejection. Similar calculations for heavier molecules indicate 
that their ejection is unlikely (Watson and Salpeter 1972). Essentially 
the coupling between vibrations and center-of-mass translational motion 
is small, because the molecular vibrational time-scale is much shorter 
than that for change in the repulsive interaction potential with the' . . 
surface (Hunter and Watson 1978). However, these calculations suffer 
from two deficiencies: they ignore transfer of energy to other decay 

channels than center-of-mass translational motion and the quantum nature 
of the energy transfer process. It should be noted that molecule forma- 
tion may occur along a pathway that is directly dissociative in the 
surface-admolecule bond (Aronowitz and Chang 1985). The general impor- 
tance of this ejection mechanism is, however, unknown, and we will 
ignore this channel in the subsequent discussion. 

The transfer of the reaction heat to other channels is very impor- 
tant. In fact, in the classical trajectory calculations referred to 
above, a somewhat arbitrary criterion had to be used to determine the 
fraction of H 2 molecules ejected, because in the absence of other decay 
channels all newly formed H 2 molecules will eventually acquire enough 
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energy in their center-of-mass translational motion to be ejected. 

Among the other energy channels that can accept the excess vibrational 
energy of a newly formed H 2 molecule are rotations of the H 2 molecule 
(rapidly followed by transfer to lattice phonons), (de)localized phonons 
of the grain, and vibrational modes of neighboring molecules. Because 
of the inherent quantum nature of the transfer process, the contribution 
of each of these processes depends on the spacing of the energy levels 
involved, as well as on the nature of the interaction. 

Except for energy transfer to neighboring molecules, the accepting 
energy channel has a much smaller energy level spacing than the initial 
vibrational energy channel. For example, a vibrational quantum of H 2 is 
3000 to 5000 cm"' (aO.4-^.67 eV), which is much larger than an H 2 rota- 
tional quantum (a 120 cm”^; a0.015 eV), a vibrational quantity of the 
f^-surface bond (a30 cm" ; a0.004 eV), or phonons (a 100 cm" ; 
a0 .01 eV). Therefore, in the energy- transfer process a large number of 
quanta have- to be created simultaneously. This process is described by 
an energy-gap law, and its rate decreases exponentially with the energy 
gap involved, that is, with the number of quanta that have to be gener- 
ated (Diestler 1974; Nitzan et al. 1975; Diestler et al. 1977). Experi- 
mental studies of tljie energy transfer from vibrational (»3000 cm” 1 ) to 
rotational (*30 cm” ) modes of small hydrides isolated in low- 
temperature matrices indeed show such an energy-gap law (Legay 1977). 

Energy transfer to neighboring molecules can be a very important 
energy decay channel for vibrationally excited molecules if the energy 
mismatch between the donor and acceptor is not too large and if the 
exchange is symmetry allowed. This energy exchange is induced by (long- 
range) multipolar (i.e., dipole or quadrupole) interactions. Theoreti- 
cally, the lifetime associated with this energy decay channel is esti- 
mated to be about 5 * 10" 12 sec for resonant conditions (Brus 1981). 
Experimental studies of resonant energy transfer of vibrationally 
excited CO indicate a lifetime of about 10" 11 sec (Legay 1977). Energy 
transfer in the Lorentzian wing of a transition is very inefficient and 
non-resonant energy transfer is, therefore, generally assisted by 
phonons (or other low-lying energy channels that match). This will 
again lead to an exponential decrease in the rate with increasing energy 
gap. This is confirmed by experimental studies (Legay 1977). Non- 
resonant energy transfer can still be very fast, if the conditions are 
favorable . 

As an example of these principles consider the different energy- 
decay channels for the excess vibrational energy of a newly formed H 2 
molecule one by one. Since this is a horaonuclear diatomic molecule, it 
does not possess a dipole moment. Near-resonant quadrupole-dipole 
transfer to neighboring molecules of the excitation in the high vibra- 
tional levels (v - 3000 cm" 1 ) therefore is probably slower by a factor 
of 10 2 to 10^ than the 10 psec time-scale for dipole-dipole interaction, 


leading to a time-scale of perhaps 10"^ to 10"® sec. It should be noted 
that interstellar grain surfaces are expected to contain many functional 
groups with vibrational modes at about 3 ym (e.g., OH and CH groups; c*\ 
§4.1), which can act as energy acceptors. Recall also that H 2 is very 
mobile on a grain surface, so that energy transfer is likely to occur 
between nearest neighbors. The mismatch in energy between donor and 
acceptor can be taken up by rotational energy of the H 2 molecule phonons 
or center-of-mass translational motion. For the lower levels 
(v - 4000 cm” 1 ) this decay channel is probably closed because of the 
large energy mismatch. 

Rotations might also directly act as energy acceptors of the excess 
vibrational energy. The experimentally measured relaxation rate of 
dipolar hydrides, which couple well with the lattice, ranges from about 
10' sec" 1 for the lighter ones (e.g., OH and NH) to 10® sec" 1 for the 
heavier ones (e.g., HC1; Legay 1977). The mass ratio of unity favors a 
high relaxation rate for this channel. However, because H 2 is nonpolar, . 
it will not couple as well to the lattice. The transfer rate of vibra- 
tion to rotation might perhaps be as fast as 10® sec" 1 . Small energy 
mismatches (i.e., less than one rotational spacing) can be taken up by 
center-of-mass translational motion of the H 2 molecule or other local 
phonon modes. Decay of the rotational energy through transfer to local 
or delocalized phonon modes (one quantum at a time) will be very rapid 
because of the similarity in energy (for example, the Debye frequency of 
H 2 0 is 133 cm" 1 ). Some of the energy may wind up momentarily in vibra- 
tional energy of the H 2 -surface bond but, again, transfer to phonons 
will be rapid and buildup of energy in this mode will not occur. Thus, 
ejection of H 2 is unlikely in this case. 

Finally, direct transfer of the vibrational energy to center-of- 
mass translational motion has been estimated theoretically to have a 
rate of about 10^ sec" 1 , but this estimate is very uncertain (Lucas and 
Ewing 1981). In this case, the excess energy in this channel is suffi- 
cient to cause ejection of a translationally "hot” (*0.5 eV; i.e., one 
vibrational quantum) H 2 molecule. 

From this discussion, it seems that for the high vibrational levels 
of H 2 , decay will occur through energy transfer to neighboring mole- 
cules. The lower levels may decay through transfer to rotational modes 
or to center-of-mass translational motion. If the latter process domi- 
nates, then ejection is likely. Clearly, laboratory studies of the 
relaxation of vibrationally excited H 2 are important to the resolution 
of this point. For heavier molecules, in particular for non-homonuclear 
molecules, energy transfer through multi-polar interaction is expected 
to dominate the reaction heat dissipation. 



5. FORMATION OF H 2 ON GRAIN SURFACES 

Ultraviolet observations have revealed that molecular hydrogen is an 
important constituent of diffuse interstellar clouds, such as those in 
front of c Oph (Morton 1975). It is generally assumed that H 2 is the 
most abundant molecule in dense molecular clouds, but direct detection 
of H 2 in quiescent clouds is difficult, because of the high UV opacity 
of such clouds and the low Einstein A values of the IR rotational tran- 
sitions. Ultraviolet-pumped, IR fluorescence from H 2 has, however, been 
observed in reflection nebulae (e.g., NGC 2023) and around HII regions 
(e.g., the Orion bar; Gatley et al. 1986; Sellgren 1986). It has also 
been observed in regions with strong molecular shocks such as the BN/KL 
region in Orion (Beckwith et al. 1978). This suggests that H 2 is indeed 
the dominant form of hydrogen in molecular clouds. 

It is difficult to form H 2 through reactions in the gas phase and, 
therefore, it is generally presumed that it is formed on the surfaces of 
interstellar grains. Because of the importance of H 2 in models of 
interstellar chemistry, the principles of grain-surface chemistry out- 
lined above will be applied in this section to the formation of H 2 in 
the interstellar medium. First, H 2 formation on bare silicate grains in 
a diffuse interstellar cloud (n Q - 10 cnT^, t - 100 K, and T d - 20 K) 
will be considered, then H 2 formation on the basal plane of graphite and 
amorphous carbon will be briefly examined. This will be followed by a 
discussion of poisoning of sites at which H atoms can be chemisorbed. 

The section will close with a discussion of H 2 formation on icy grain 
mantles inside molecular clouds. Hydrogen formation on organic refrac- 
tory dust is presumably very similar to this latter case. 


5.1 


H 2 Formation on Silicate Surfaces 


Hydrogen atoms on a silicate grain can be chemisorbed as well as physi- 
cally adsorbed on top of a chemisorbed layer (cf. §4.2). There are 
about 10° of each of these sites on a grain with a radius of 10^ A. 

Under the conditions in the interstellar medium, migrating physically 
adsorbed atoms will react with chemisorbed atoms forming H 2 , which will 
evaporate quickly. With a sticking coefficient of unity and the gas 
parameters stated above, the H accretion time-scale, t a , on a 1000 A 
grain is about 500 sec. Thermal evaporation from chemisorption sites is 
completely unimportant. For physical adsorption sites the evaporation 
time-scale is about 2 x 10 1 * sec. The migration time-scale for physi- 
cally adsorbed atoms is, however, only 3 * 10 "^ sec and, if there is an 
empty chemisorption site, the physically adsorbed atom will find it 
before it evaporates. Most H atoms will, thus, be bound in chemisorp- 
tion sites. Occasionally, a chemisorbed atom can tunnel to an adjacent 
site on a time-scale T t of about 10^ sec, assuming a barrier with a 
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height of 1 eV and a width of 1 A. If this site is empty, the H atom 
will be chemisorbed again until it tunnels on. If there is another 
H atom present inthis site, the entering H atom will only be physically 
adsorbed. It may react with the H atpm already present or migrate on. 
With an activation barrier of 1000 K and a width of 1 A, the time-scale 
for reaction, r p , is 1.3 x 10”° sec. This should be compared to the 
time-scale for migration (3 * 10"^ sec). The probability of reaction, 
p r , is therefore only 2 x 10"^ and the H atom will migrate about 500 
sites before it reacts (if it is not trapped in an empty chemisorbed 
site before that). The probability for evaporation during this migra- 
tion and reaction time is given by T p /x ev ; for the parameters above, 
this is negligible, and all accreting H atoms will form Hj. 

When the dust temperature is increased, evaporation of H atoms in 
the physically adsorbed layer will become important. A critical temper- 
ature, T cp , can be defined for which the probability for reaction with a 
chemisorbed atom is equal to the evaporation rate. This critical tem- 
perature is given by 


T = E. /mix v ) 
cr b r z 


( 16 ) 


where is the binding energy and v z the characteristic vibration 

frequency of the H atom in the physically adsorbed layer (cf. 
eq. [6]). For the values given above, this critical temperature is 
50 K. Note that this critical temperature is somewhat lower than that 
calculated by Hollenbach and Salpeter (1971). They envisioned a situa- 
tion with a (small) number of enhanced binding sites interspersed in an 
otherwise physically adsorbed layer. However, they neglected the acti- 
vation barrier for reaction. In the analysis above, the activation 
barrier for reaction is crucial. If the activation .barrier for reaction 
were actually less than the binding energy of the H atom in the physi- 
cally absorbed site, then H 2 formation would occur until evaporation 
from the chemisorbed layer dominated the loss of H atoms. The labora- 
tory studies measure, however, both the physical adsorption energy and 
the activation energy rather directly, and the latter is indeed larger 
than the former. Above the critical temperature, the reaction rate will 
drop by a factor T ev /T p 3 3 x 10”' exp(750/T^) ; at 75 K, this factor is 
6 x 10”^. Finally when the dust temperature is very high and evapora- 
tion of physically adsorbed H atoms is very rapid, only those impinging 
gas-phase atoms that have enough kinetic energy to directly overcome 
the activation barrier for reaction with the chemisorbed atoms will 
react.' The rate of H 2 formation is then reduced by a factor 

exp(-E a /kT) * 0.25 exp(1000/T) a 10 . Here the fraction of 

a site covered by an H atom, takes into account the smaller cross sec- 
tion of an H atom compared to the site area. Note that evaporation of 



chemisorbed H atoms is unimportant until the grain temperature reaches 
about 800 K. 

5.2 H 2 Formation on Graphite Surfaces 

The formation of H 2 on the basal plane of graphite in the diffuse inter- 
stellar medium is very similar to that on silicates. As on silicates, 
atomic hydrogen can be both chemisorbed and physically adsorbed on 
graphite with binding energies of about 2 eV and 0.07 eV, respectively 
(see §4.2). The most important parameter for H 2 formation, the activa- 
tion barrier for the reaction between physically adsorbed and chemi- 
sorbed H, is also very similar for the two surfaces (*1200 K; King and 
Wise 1963; versus =*1000 K; Wood and Wise 1962). The only difference is 
the presence of so-called "active surface areas" on graphite (i.e. , the 
edges of the basal plane). Hydrogen atoms can attack the carbon skele- 
ton at those areas and form saturated hydrocarbons such as methane. 

There is a considerable spread in the activation barriers reported for 
this set of reactions (1700-5000 K; King and Wise 1963; Wood and Wise 
1969; Rye 1977; Balooch and Olander 1975), which may be due to differ- 
ences in the surface structure of the samples (Veprek and Haque 1975). 

In any case, the activation energy is larger than that for reaction of a 
migrating H atom with a chemisorbed H atom on the basal plane; as a 
result, diffusing H atoms generally form H 2 on the basal plane instead 
of reacting on the edges. This is even more pronounced in the inter- 
stellar medium because H coverage on the basal plane will be high. In 
this respect, it should be noted that in the one series of experiments 
in which the H atom coverage on the basal plane was high, H 2 was the 
dominant product at low temperatures (King and Wise 1963). Draine 
(1979) has argued that a migrating H atom would abstract a chemisorbed 
H atom (part of a CH^ group) on these edge structures to form H 2 rather 
than form CH^ at low temperatures. However, this H abstraction reaction 
is likely to have considerable activation energy (*8000 K; Balooch and 
Olander 1975), and from the experiments we infer that it does not play a 
role until the surface temperature is much higher ( *800 K). At that 
temperature, CH^ becomes thermodynamically unstable, and the equilibrium 
will shift toward hydrogen and (solid) carbon (Wood and Wise 1969). 

5.3 Poisoning of Chemisorbed Sites 

We will now consider poisoning of the chemisorbed sites by other 
species. That this is important can be seen by comparing the accretion 
rate with the thermal evaporation rate. For small grains (<100 A) or 
warm grains (T d > 22 K), physically adsorbed H atoms evaporate before 
the next H atom accretes, and H 2 formation through recombination of two 
physically adsorbed H atoms is inefficient. In the interstellar medium, 



poisoning will occur when another atom blocks the access of H to the 
chemisorbed site. The poisoner may be physically adsorbed as long as it 
is essentially immobile. Note that if there are still a few chemisorbed 
sites available on a grain that H 2 formation will occur by the mechanism 
outlined above. Poisoning is, thus, only effective if almost all of the 
sites are covered by poisoning atoms. The time-scale for poisoning is 
then the accretion time-scale per site of poisoners, which is about 
5 x 10^{10 -1 */X) yr, with X the relative gas-phase abundance of the 
poisoning atoms. For oxygen, nitrogen, and carbon atoms this time-scale 
is very short compared to the lifetime of the cloud or the H 2 photodis- 
sociation lifetime, and they might therefore be very effective poison- 
ers; this is even more likely because they are also chemisorbed to these 
sites. However, H atoms will be able to tunnel through the activation 
barrier for reaction and form H 2 O, NH^, and CH^, which is energetically 
favorable. However, in the diffuse interstellar medium, photodesorption 
is an effective way of cleaning a surface of these molecules, since 
their photodesorption time-scale in an unshielded region is only about 
700 yr (see §9.1). 

Poisoning by other atoms might also be important in the diffuse 
medium. The observed total surface area of interstellar dust is about 

3 x ICT^ 1 cm^ per H atom. This number might be somewhat larger (50%) if 

the grain size distribution extends to below 100 A. With a surface area 
per site of about 7 , this dust surface area corresponds to about 

4 x 10“® sites per H atom. Thus, species with a total gas-phase abun- 
dance in the diffuse interstellar medium comparable to this number could 
poison the chemisorbed sites on a grain surface. Note that we have to 
correct for the underlying depletion caused by the incorporation of some 
elements in refractory dust grains produced by stars. This implies that 
poisoning by atoms such as S, Si, Na, Mg, Fe, Al, Ca, and Ni can be of 
importance. Some of these (e.g., S and Si) will form volatile hydrides 
which are probably also easily photodesorbed, whereas others (i.e., Na 
and'K) may be ejected upon electron recombination (cf. §6). The transi- 
tion metals, as well as Ca, will form hydroxide groups on interstellar 
grain surfaces (§6). Presumably, H will also be chemisorbed on such 
surfaces. Moreover, although their abundance might be sufficient to 
completely poison the chemisorbed sites on a time-scale of only about 
10 7 yr, this assumes that these species are used as efficiently as 
possible (i.e., one atom per chemisorbed site). This is unlikely, since 
these species are immobile at 20 K. It seems, therefore, that poisoning 
of grain surfaces cannot stop H 2 formation in the diffuse interstellar 
medium. 



5.4 H 2 Formation on Grain Mantles 

Finally, H 2 formation on icy surfaces will be considered. Inside dense 
molecular clouds, interstellar grains. will accrete mantles consisting of 
such simple molecules as H 2 0, NH3, CO, and CH3OH (see §7). Such icy 
grain mantles have only physical adsorption sites, with a typical bind- 
ing energy of about 350 K; this is owing to the high surface coverage of 
H 2 (Govers et al. 1980; Helens and Hagen 1982; see also §4.2). The 
thermal evaporation rate of H is then about 2 * 10 2 /sec, which is lower 
than the accretion time-scale of H atoms on a grain (-10^ sec for a 
1000 A grain and a typical H density of 1 cm" 3 ) and H 2 formation by 
direct recombination of physically adsorbed H atoms on a grain surface 
will not occur. Of course, by the time a grain mantle has been formed, 
hydrogen is already completely converted into the molecular form through 
reactions on bare graphite or silicate grain surfaces. 

Although icy grain mantles do not have chemisorbed sites, molecules 
such as H 2 C0, H 2 S, N 2 H 2 , N 2 H^, and metal hydrides can act as enhanced 
binding sites for H 2 formation (Tielens and Hagen 1982). During its 
lifetime on a grain surface, an H atom may tunnel through the activation 
barrier associated with H abstraction from such molecules producing H 2 
and leaving a reaction site <cf. §4.5). Thus, H 2 can also be formed and 
reformed inside dense molecular clouds'. It should be stressed that 
because radical -radical reactions are activationless, if there are 
radicals on the surface, then an H atom will react with them rather than 
tunnel through an activation barrier and form H 2 . Reaction of H with 
non-radicals such as 0 2 and CO can also be important in this respect. 
Thus, when the gas -phase atomic H abundance decreases below the abun- 
dance of heavier species such as 0, CO, and 0 2 , H 2 formation on grain 
surfaces will stop (Tielens 1983b). Typically, this occurs at a density 
of about 3 * 10 4 cm -3 (inert conditions; see §7). This will, however, 
not increase the abundance of H in the gas phase, because the H atom 
will now be incorporated in such molecules as H 2 0, NH^, and H 2 C0. 


6. DEPLETION OF THE ELEMENTS 

UV absorption line studies have shown that the gas-phase abundance of 
many elements in the diffuse interstellar medium is less than expected 
from cosmic abundances. The degree of depletion varies from element to 
element, as well as from one line of sight to another (see chapter by 
Jenkins, elsewhere in this volume). The missing elemental fraction is 
generally thought to be locked up in grains and, thus, elemental abun- 
dance studies may yield information on the composition and physical and 
chemical properties of interstellar grains. Several processes can 
influence the gas-phase abundance of an element. Clearly, dust 
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condensation in the high-temperature and high-density outflow from 
stellar sources (Field 1974) is important, as evidenced, for example, by 
the large fraction of silicon locked up in silicates (cf. §1.1). 
Destruction of dust in interstellar shocks is also important, given the 
relatively high gas-phase abundance of elements such as Ca, Fe, and Si 
in high-velocity clouds (i.e., recently shocked clouds) compared to the 
general diffuse interstellar medium (Spitzer 1976). Finally, accretion 
of atoms in the interstellar medium is also an important process (Snow 
1975), since elemental abundances lower than 10 -1 times the cosmic 
abundance are difficult to achieve by stardust formation alone (see 
chapter by Jura, elsewhere in this volume). The complexity of these 
interplaying processes has defied a general theory of elemental deple- 
tions in the interstellar medium. Considering the scope of this book, 
in this section we merely want to point out some considerations concern- 
ing the third process, the accretion of elements in the diffuse inter- 
stellar medium. 

Nitrogen, oxygen, and sulfur are expected to be little depleted in 
the interstellar medium, because their main condensation products around 
stars are very volatile (e.g., NH^, H 2 0, and H 2 S). Although these 
elements are presumably chemisorbed on interstellar graphite and sili- 
cate surfaces, they are expected to react in the diffuse medium with 
atomic hydrogen on grain surfaces and thereby form their volatile 
hydrides. These are probably easily photodesorbed (cf. §9). Once they 
are in the gas phase, photodissociation will rapidly destroy these 
molecules. Sputtering in low velocity («20 km/sec) shocks can also be 
of importance for the removal of these saturated molecules. Some oxygen 
will be locked up in refractory silicate grains and, possibly, in 
organic refractory material. 

Carbon takes a special place among the more abundant heavy ele- 
ments. Not only is solid carbon the thermodynamically favored phase at 
high temperatures and densities (i.e., in the outflow froms stars), but 
carbon can also form large molecules that can be fairly resistant to 
destruction. Observations indicate that a considerable fraction of the 
carbon can be locked up in grains made up of cross-linked macromolecules 
(cf. §1.3). The physical processes important in the formation of the 
organic refractory dust component in the interstellar medium are dis- 
cussed in §8. Here we merely emphasize that carbon can be highly 
depleted in the interstellar medium. 

Silicon and phosphorous can form a refractory high-temperature 
condensate in the outflow around stellar objects (i.e., MgSiO^ and Fe^P; 
Fegley and Lewis 1980). In the diffuse interstellar medium, it is 
expected, however, that both atoms when accreted on a grain surface will 
form saturated, volatile hydrides (i.e., SiH^ and PH^) in a manner simi- 
lar to that of carbon and nitrogen. Again, photodesorption will prob- 
ably prevent large depletions in the diffuse medium. Silicon can 
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substitute for carbon in organic molecules, and its depletion in the 
interstellar medium may, therefore, vary through incorporation into the 
organic refractory dust component formed in molecular clouds (cf. §8). 

The depletion pattern of the alkali and alkaline metals is somewhat 
puzzling in view of their closeness in the periodical system. Sodium, 
potassium, and magnesium are little (about a factor of two) depleted in 
stardust (presumably composed of silicates) and do not show large dif- 
ferential depletion effects in the diffuse interstellar medium. In 
contrast, Ca is highly depleted and shows large variations. Both Na and 
K have ionization potentials below the grain work function and, thus, 
charge effects will not complicate their sticking (i.e., S » 1; cf. 

§4.3). The positive ions Na* and K + have an inert-ga3 electronic 
structure, which inhibits their reaction with H upon accretion. How- 
ever, once accreted an electron is expected to be able to tunnel through 
the work function and recombine with the alkali metal ion. If this 
recombination occurs in a Rydberg state (high principal quantum number), 
the extended electron cloud coupled with the repulsive interaction with 
the surface will then, presumably, lead to desorption (see the discus- 
sion on photodesorption in §9.1). Recombination may also occur directly 
in an antibonding valence orbital of the surface-adion bond. If this 
scenario is correct, then one would expect that because of their similar 
electronic structures and low ionization potentials, Li, Rb, and Cs are 
also little depleted in the interstellar medium. 

Because the ionization potentials of Mg and Ca lie above the work 
function of interstellar grains, sticking might be complicated by recom- 
bination effects (cf. §4.3). The large depletion of Ca suggests, how- 
ever, that electrons in interstellar grains are not very mobile on the 
short collision time-scale (*10 -1 3 sec ). Since the cations do not have 
an inert-gas electron structure, they probably react immediately upon 
accretion with atomic hydrogen on the surface, which may actually pre- 
vent their ejection upon electron recombination. Both elements will 
react with hydrogen to form the dihydrides. However, -whereas MgH 2 has 'a . 
covalent structure, CaH 2 is a salt (Mueller, Blackledge, and Libowitz 
1968). Presumably, MgHg is therefore easily photodesorbed and CaH 2 is 
not. If this is indeed the essential difference, then it is expected 
that other elements that form covalent hydrides, such as Be, Zn, and Cd, 
would be little depleted (Zn is discussed below). At this point, it 
should be recalled that hydrogen abstraction reactions are very impor- 
tant on interstellar grain surfaces. The Ca hydride remaining on the 
surface will therefore react often with H (i.e., CaH 2 -*■ H CaH + H 2 , 
followed by CaH + H ♦ CaH 2 ). This cycle will continue until CaH is 
poisoned by another atom, presumably 0, leading to hydroxide forma- 
tion. It is expected that the Ca(0H) 2 will be bonded quite strongly to 
the surface. 
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The transition metals Ti, Mn, Fe, and Ni are highly depleted in the 
interstellar medium. In contrast, the depletion of Zn is quite small. 
The special feature of the electronic structure of the transition metals 
that presumably lead to this depletion pattern is the presence of 
unfilled d orbitals, which can form strong complexes with electron-rich 
material (Pauling 1964). The unfilled 3d orbitals available should 
provide strong bonds to a grain surface. These elements can also react 
with atomic H on the surface to form mono- and dihydrides, which could 
be stabilized by the formation of a coordination complex with a ligand 
group on the surface (e.g., highly polar or polarizable groups with 
unshared electron pairs such as an OH group in a si land moiety on a 
silicate surface or the ir electron system in PAHs, amorphous carbon 
particles, or graphite). This complex will still be strongly bound to a 
grain surface and photodesorption will be of little importance. In this 
respect it is important to note that Zn has all its 3d orbitals filled 
and does not readily form coordination complexes. Zinc hydride has a 
covalent character (Mueller, Blackledge, and Libowitz 1968), and is 
probably easily photodesorbed or sputtered. Cadmium, which has an 
electronic structure similar to that of Zn, also forms a covalent 
hydride and should, therefore, have a similar depletion pattern. 

In summary, the differential elemental depletion patterns in the 
diffuse interstellar medium can be ascribed to the chemical properties 
of the elements. Of particular importance is the ability to form 
hydrides and coordination complexes. Photodesorption and electron 
recombination may play a key role in preventing depletion of some ele- 
ments. In contrast to Barlow (1978), we do not think that ejection upon 
reaction is an important process (cf. §4.6). We also caution against 
the straightforward use of surface bond energies of the elements based 
on chemisorption theory to explain observed depletion patterns. Chemi- 
sorption is a very complex process and is far from theoretically under- 
stood (Einstein et al. 1980). 


7. COMPOSITION OF ICY GRAIN MANTLES 

The composition of interstellar grain mantles has been calculated, 
employing gas-phase as well as grain-surface reactions (Tielens and 
Hagen 1982; Tielens 1983b; d'Hendecourt, Allamandola, and Greenberg 
1985). Figure 13 shows the composition of the gas phase and the grain 
mantle as a function of the gas-phase hydrogen density (Tielens and 
Hagen 1982). In the subsequent discussion we will concentrate on the 
global composition of the grain mantles and ignore to a large extent the 
abundances of trace species. 

First, we will concentrate on the gas-phase composition of molecu- 
lar clouds (Fig. 13a), since it sets the stage for the diffusion-limited 
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Figure 13. Molecular composition of interstellar clouds as a function 
of density in the gas phase (a) and grain mantle (b). 


chemistry on grain surfaces. Hydrogen is mainly in molecular form 
(H2). The atomic hydrogen abundance is inversely proportional to the 
total density', because the H2 formation time-scale depends on the 
(grain) density, whereas the destruction time-scale (by cosmic rays) 
does not. The most abundant molecule in the gas phase after H 2 is CO, 
which contains most of the elemental carbon in the gas phase. Oxygen 
and nitrogen are either in atomic or molecular form (0 2 and N 2 ). Note 
that H 2 0 contains only a small fraction of the gas-phase oxygen. These 
results are very typical for time- independent gas-phase chemistry 
calculations or for time-dependent gas-phase chemistry calculations at 
late stages (>10° yr; Mitchell et al. 1978; Prasad and Huntress 1980; 
Graedel et al. 1982). 

The calculated grain mantles consist mainly of the molecules H2O, 
H2CO, Nj , O2, CO, CO2, HgOj, and NH3 and their deuterated counterparts 
in varying ratios. The exact composition depends strongly on the condi- 
tions in the gas phase. Depending on the gas-phase composition, three 
different physical regimes can be discerned, a reducing, an inert, and 
an oxiding "atmosphere." For low densities (n < 5 * 10^ cm"^), atomic 
hydrogen is more abundant in the gas phase than the heavier species 
(e.g., 0, O2, and CO) and the gas is largely reducing. As a result, 
most species will be hydrogenated on the grain surface (e.g., H2O and 
H2CO). Unsaturated molecules, such as CO, are trace species in the 
grain mantle, whose abundance depends directly on their gas-phase abun- 
dance relative to that of atomic H and is, thus, a strong function of 
the total hydrogen density. An exception should be made for N2, which 
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does not react with atomic H. Ammonia is therefore a trace species, 
whose abundance depends on the gas-phase abundance of atomic N. 

At higher densities (n > 5 * cm”3), the atomic H density in the 

gas phase is less than that of the heavier species and most accreting 
species will not be hydrogenated. The resulting grain mantle depends 
then on the 0 to O2 ratio in the gas phase, since 0 is the next most 
abundant, mobile radical. If this ratio is low, then only few grain 
surface reactions will occur (i.e., inert atmosphere). The grain mantle 
will then mainly reflect the gas-phase composition (e.g., CO, Oj, and 
Ng). If the atomic 0 to O2 ratio in the gas phase is high (oxidizing 
atmosphere) then atomic 0 and CO will react; C0 2 will then be the domi- 
nant grain-mantle constituent. In both cases hydrogenated molecules 
such as H 2 0 and H 2 CO are only trace species, whose abundance depends on 
their abundance in the gas phase and on the ratio of atomic H to heavier 
species. 

Infrared observations show that interstellar grain mantles consist 
mainly of H2O and CHgOH; CO and NH3 are only trace species (Tielens 
et al. 1986a). This implies that grain mantles are formed in a "reduc- 
ing" atmosphere and, thus, at relatively low densities (»10^ cm"3). 

Note that CH^OH may result from H addition reactions to F^CO. There is 
some experimental justification for the occurrence of this reaction at 
low temperatures (Tielens and Hagen 1982). It has, however, not been 
included in the calculations described above. Finally, in some cases 
evidence for grain mantle formation in a more inert atmosphere has also 
been observed (Tielens et al. 1986a). 

It is interesting to consider the deuteration of species in a grain 
mantle. Figure 14 shows the relative abundance of a deuterated species 
with respect to its hydrogenated counterpart (Tielens 1983b). Gener- 
ally, this ratio reflects the ratio of atomic deuterium to atomic hydro- 
gen in the gas phase, essentially because the grain-surface reaction 
networks are very similar. Thus, this ratio increases strongly with 
increasing gas-phase density owing to the decrease in the gas-phase 
abundance of atomic hydrogen (cf. Fig. 13a). It is important to note 
that although most of the hydrogen is in molecular form, most of the 
deuterium is actually in atomic form in the gas phase, and the ratio of 
atomic deuterium to atomic hydrogen is much larger than expected from 
elemental abundances (Tielens 1983b). Molecular hydrogen is destroyed 
in the gas phase by cosmic ray ionization leading eventually to H^. 
However, in its destruction network, HD competes with other molecules 
such as CO and N2 for reaction with H^ forming H2D'*\ Reaction of this 
molecular ion with neutral molecules such as CO and N2 (leading to the 
deuterated ions DCO + and N2D + ) , followed by dissociative electron 
recombination produces atomic deuterium. On the basis of their relative 
abundances, one out of about every 20 cosmic ray ionizations of H2 will 
lead to HD destruction; thus, the HD destruction rate is comparatively 
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Figure 14. Concentration of the most abundant deuterated molecules 
relative to their hydrogenated counterparts in the grain mantle as a 
function of the total hydrogen density in the gas phase. 

much larger than that of H 2 (Tielens 1983b). In fact the high deu- 
teration observed for molecular ions in molecular clouds implies that 
the D/H ratio is much larger than the elemental abundance ratio. Note 
that the grain-mantle abundances of some molecules (e.g., HDCO, D 2 C0 and 
to a lesser extent HDO) do not directly reflect the atomic D/H ratio in 
the gas phase. These result from differences in the grain-surface reac- 
tion network for these species (e.g., H abstraction reactions are of 
importance; Tielens 1983b). 

8. UV PHOTOLYSIS, TRANSIENT HEATING, AND ORGANIC REFRACTORY DUST 

Greenberg and his co-workers and others (Greenberg 1963; Greenberg and 
Yencha 1973; Hagen et al. 1979; Greenberg 1979; d'Hendecourt et al. 

1986; Donn I960; Donn and Jackson 1970; Khare and Sagan 1973) have 
suggested that UV photolysis of interstellar grain mantles inside molec- 
ular clouds can produce complex molecules. These complex molecules are 
not as easily destroyed by interstellar shock waves as H 2 0 ice (Draine 
and Salpeter 1979) and could conceivably form an important component of 
the dust in the diffuse interstellar medium. Processing by cosmic rays 
might play a similar role (Moore and Donn 1982; Strazulla et al. 

1983). Infrared spectroscopy of protostars deeply embedded inside dense 
molecular clouds suggests that the composition of interstellar grain 
mantles is more complex than can be explained by grain-surface reactions 
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alone (Lacy et al. 1984; Helens et al. 1984, 1986a). This may indicate 
the importance of energetic processing by UV photons or cosmic rays. 

It should be stressed that grain-surface reactions alone only for:.i 
simple, saturated hydrides (cf. §7). The formation of the organic 
refractory dust component, consisting of long carbon chains (i.e., 
polymeric carbon), in such a reducing atmosphere as the interstellar 
medium is difficult to understand. This section will discuss this 
problem and will identify the two key processes — that is, UV photolysis 
producing radicals and radical diffusion triggered by transient 
heating — that may convert interstellar grain mantles consisting of 
simple saturated hydrides into the complex organic refractory dust 
component of the interstellar medium. 

In contrast to grain-surface reactions occurring during accretion 
of gas-phase species, UV photolysis mainly takes place in the bulk. 
Ultraviolet photolysis of interstellar grain mantles can produce reac- 
tive- radicals, often in electronically or vibrationally excited states 
(Hagen 1982). These "hot" radicals can react with their neighbors or 
they can diffuse through the matrix and react with other species present 
(Hagen, Allamandola, and Greenberg 1979). Because of the large amount 
of excess energy available, even reactions with molecules having closed 
electron shells, which normally have appreciable activation barriers, 
might take place. Because, in general, reactions only take place when 
the reactants approach along well-defined reaction paths and because 
collisions of the hot reactant and the host material will transfer 
energy to the matrix (e.g., heat), reactive species can also be stored 
at low temperatures in the grain mantle. Upon warm-up those radicals 
can start diffusing through the host material and react with each 
other. Thus, during photolysis and warm-up complex molecules can build 
up, producing an organic refractory mantle on interstellar dust grains 
(Greenberg 1979). 

The chemical processes taking place during grain mantle photolysis 
are complex. Often there are several photo-decomposition channels of a 
primary molecule. The photoproducts can also react with many of the 
neighboring unsaturated molecules. Moreover, the resulting secondary 
molecules themselves can often also be photolyzed, leading to a third 
generation of molecules. This complexity is compounded by the loss of 
volatile components owing to temperature spikes generated by grain-grain 
collisions (d'Hendecourt et al. 1982) or cosmic rays (Leger et al. 

1985). The complexity of these processes has inhibited theoretical 
studies of photolyzed grain mantles. Extensive and painstaking labora- 
tory studies are required to understand these processes. Although 
some laboratory studies have already been undertaken (Hagen 1982; 
d'Hendecourt et al. 1986; van de Bult 1986), this field is still in its 
infancy. Little more beyond the general outline of the chemical scheme 
(UV photons produce radicals which form complex molecules) can presently 
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be said. Among the typical radicals produced by the photolysis of H 2 0- 
rich ices are HCO, HOCO, and OH. Simple molecules that have been pro- 
duced by photolysis of such mixtures are HgCO and C0 2 (Hagen 1982; 
d'Hendecourt et al. 1986). A general conclusion seems to be that ali- 
phatic rather than aromatic hydrocarbon mixtures are readily formed. 
Ultraviolet photolysis of oxygen-rich mixtures generally leads to the 
formation of residues containing large amounts of carboxylic acids. If 
NH^ is also present in the initial mixture, then amino groups are also 
formed (Hagen, Allamandola, and Greenberg 1979; d'Hendecourt et al. 
1986). 

Let us consider quantitatively the conversion of interstellar 
grain mantles into organic refractory dust by means of UV photolysis. 


The absorption cross section of H 2 0 ice below 1600 A is about 
10^ 1 mole -1 cm -1 , about 4 times larger than for gas-phase H 2 0. In an 
unattenuated interstellar FUV radiation field (10° photons cm -2 sec -1 ; 
Habing 1968), the number of photons absorbed by one H 2 0 molecule in a 
grain mantle, N ab 5 < # v = 0), is then about 2 x 10 - ^ sec -1 . Taking a UV 
dust albedo of 0.6, a mean cosine of the scattering angle of the dust of 
0.7 (Lillie and Witt 1976), and a UV dust absorption cross section 
3 times the visual extinction, A y , the number of UV photons absorbed by 
one H 2 0 molecule in a grain mantle as a function of depth in a semi- 
infinite cloud is given by (Flannery et al. 1980) as 


N abs ( V 3 2 * 10 ‘ 9 e*P(-2A v > 


(17) 


Thus, deep inside dense clouds (A y > 7 mag) UV photolysis by the 
ambient interstellar radiation field is unimportant over the lifetime of 
a molecular cloud (3 x 10^ yr). The UV photon field inside a dense 
cloud may actually be larger than estimated from this treatment because 
of the presence of local UV sources. For example, cosmic ray excitation 
of H 2 may lead to a mean intensity of the UV radiation field inside the 
cloud of 10^ photons cm -2 sec -1 (Prasad and Tarafdar 1983). Thus, for a 
visual extinction larger than about 6 mag this field will dominate. 

There may also be other sources of UV photons inside dense clouds, such 
as shocks or newly formed stars, whose flux is difficult to estimate but 
which might be as high as 5 * 10 J photons cm" 6 sec" (Norman and Silk 
1980). Taking 10^ photons cm -2 sec" 1 as the value for the mean inten- 
sity of the UV radiation field inside dense clouds yields, for the 
number of UV photons absorbed by one H 2 0 molecule in a grain mantle, 

N a bs = 2 x 10 ~ 14 sec" 1 . Thus, over the lifetime of a molecular cloud 
every H 2 0 molecule in a grain mantle will have absorbed about 20 UV 
photons . 

Not all of the absorbed UV photons will, however, lead to radical 
production. The efficiency, f uv , of the photodissociation reaction 
(number of radicals produced per UV photon absorbed) depends on the 
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particular electronic state involved. In the gas phase, photodissocia- 
tion is the dominant channel (Okabe 1978; Calvert and Pitts 1966); in a 
solid, many other decay channels compete. These include photo-assisted 
reactions with neighbors, and dissipation of the energy in the phonon 
modes of the solid (i.e., heat; Bass and Broida 1966; Bondybey and Brus 
1980). The accessibility of the various channels depends on the compo- 
sition of the ice, the identity of the molecule that has absorbed the UV 
photon, and the particular state it is excited into. Considering HgO 
under a few well-defined conditions, several pathways can be identified 
and followed. The dissociated H 2 0 molecule may recombine immediately; 
the "hot" product may recombine with a neighbor to form a new species 
and it may also recombine with another radical to reform the original 
molecule. In that case, the UV photon has merely promoted diffusion of 
a previously stored radical. The importance of this latter process 
depends on the average distance between stored radicals as compared to 
the average distance a hot product will diffuse before it is trapped. 
Note that this process will limit the concentration of radicals that can 
be stored in a grain mantle. Laboratory studies of UV photolysis of H 2 0 
and CO matrices yield values for f ranging from 0.25 to 0.005 when 
the H 2 0/C0 ratio increases from 10"" to 2 (Hagen 1982). In the former 
case, the H and OH produced will react with a CO, producing HC0 and 
H0C0. In the latter case, there is a much greater probability of reac- 
tion of the H with an OH previously produced, thereby reforming H 2 0. 

The cage effect might limit f uv , thus, to values less than 0.25. For 
the lower CO concentrations, the efficiency is limited by the recombina- 
tion of the hot product reforming HgO. 

The limiting concentration of H 2 0 photodissociation products in a 
pure H 2 0 matrix is estimated to be about from these experiments 
(Hagen 1982). In photolysis experiments involving more complex mixtures 
designed to duplicate an interstellar grain mantle, the conversion 
efficiency also falls in the 1% range (d'Handecourt et al. 1986). Water 
is a major component of interstellar grain mantles (Tielens and Hagen 
1982; Tielens et al. 1986a). On the basis of these experiments, the UV 
photolysis efficiency for interstellar grain mantles is about 0 . 5% - 
With this efficiency, about 10 % of the mantle molecules will have been 
converted into photolyzed products over the lifetime of the cloud. 

Thus, while photolysis will not change the composition of accreted grain 
mantles drastically it can increase the complexity of the molecules in a 
grain mantle. It should be emphasized, however, that most of these 
photolyzed products will not be very large. In particular, photolysis 
of H 2 0 will lead to reactions of its photolysis products, H, 0, and OH 
with neighboring simple molecules such as CO, thereby forming C0 2 , H 2 C0, 
and CH^OH. In an interstellar grain mantle, molecules containing CH 2 , 
CH^, and NH 2 groups can also be photolyzed, adding to the complexity of 
the mixture. However, most product molecules are expected to be 
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volatile and will not survive long in the harsh environment of the 
diffuse interstellar medium. The efficiency of the formation of very 
complex molecules is expected to be low. Using the values quoted above > 
for a molecule containing N building blocks, f « 20(0.005) n , where 
n = 2 log N. UV photolysis alone can, therefore, not account for the 
presence of an abundant complex organic refractory dust component in the 
interstellar medium. 

Laboratory studies show that after warm-up of a photolyzed icy 
mixture, a complex molecular mixture remains. Presumably, this mixture 
has been partly formed by radical diffusion and reaction during the 
warm-up (the exploding grain mantle; see §9.1). However, the efficiency 
of this process has not been estimated. For the sake of argument we 
will assume that It of all the radicals present at the warm-up will form 
such complex molecules (i.e., 10“ 4 of all the molecules in the grain 
mantle). The equilibrium concentration of radicals in an ^O-rich 
matrix (*1%) is attained on a time-scale of about 10® yr inside a dense 
cloud, which is fast compared to the accretion time-scale (10^/n Q yr). 
Warm-up of the grain mantles may be due to grain-grain collisions 
(Greenberg 1979) with a time-scale of about 10°/(n Q ) 1/2 yr in molecular 
clouds (Draine 1984). It has also been suggested that warm-up results 
from transient heating by cosmic rays with a time-scale of about 
3 x 10 4 yr (see §9.1). The limiting step seems, thus, to be the accre- 
tion on grains. Taking an accretion time-scale of 10° yr, appropriate 
for a density of 10® cm“® and a cloud lifetime of 3 * 10^ yr, implies 
that each atom cycles about 30 times through a grain mantle and thus has 
0.3% probability of being incorporated into the organic refractory 
dust. Since an atom cycles about 50 times through a molecular cloud 
during its lifetime in the interstellar medium, there is a 15% probabil- 
ity that it will be part of the organic refractory dust component during 
some part of its lifetime. 

It should be emphasized that there are two key steps in models of 
the UV photolysis of icy grain mantles and the production of complex 
organic refractory material. First, such models hinge on the presence 
of sources of UV radiation inside a dense cloud. Given a cloud lifetime 
of 3 x 10' yr, the ambient interstellar UV radiation field is unimpor- 
tant for UV photolysis when A v is larger than about 7 mag. Second, 
transient heating of grain mantles is required to bring the stored 
radicals together for reaction. The efficiency of such explosions is, 
however, unknown. A conversion efficiency of about 1% of the stored 
radicals is required to explain the presence of about 25% of the carbon 
in the form of organic refractory dust. 
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9. GRAINS AND GAS-PHASE COMPOSITION OF MOLECULAR CLOUDS 

One of the most pressing problems remaining in theoretical studies of 
the gas-phase composition of molecular clouds is the contribution of the 
grain chemistry. Grains can act as sources and as sinks for gas-phase 
species. The importance of grains can be easily illustrated by recall- 


1978). Thus, for n Q = lO^cm"^ gas-phase molecules accrete much 
faster on a grain than a cloud collapses to form a star. One might 
therefore expect that no molecules will be observed at such high densi- 
ties in molecular clouds (Watson 1976; Greenberg 1979). This is in 
stark contrast, however, with the observations. Molecules such as H 2 CO 
and NHg are routinely observed in dense cloud cores. Moreover, elemen- 
tal depletions in grain mantles range up to only 40% and are typically 
much less (Tielens et al. 1986a). This suggests that there is an effi- 
cient exchange between the gas phase and grain mantles. However, 
although the accretion time-scale may presently be less than the free- 
fall time-scale, it need not have been in the past. In particular, a 
newly formed star will raise the temperature of the dust in its neigh- 
borhood above the critical temperature for grain-mantle formation 
(>30 K, compare the accretion time-scale of abundant species with their 
evaporation time-scale; see Table II). 

The high abundance of molecular hydrogen in molecular clouds also 
suggests an important contribution from grain-surface chemistry to the 
gas-phase composition. There is no efficient gas-phase reaction network 
for the formation of Hj. For that reason molecular hydrogen is gener- 
ally assumed to be formed on grain surfaces by recombination of two 
hydrogen atoms or by hydrogen abstraction reactions (see §5). This 
assessment is, however, based on negation of other possibilities and 
there is no direct (chemical) evidence. Although the formation and 
ejection of H 2 on grain surfaces does not necessarily imply that other 
grain surface reactions contribute to the gas-phase composition, it is 
suggestive. Certainly, there is good observational evidence for the 
formation of molecules by grain-surface reactions in the composition of 
grain mantles (Tielens and Hagen 1982; Tielens et al. 1984; Tielens 
et al. 1986a). 

As argued in §8, grain-mantle "explosions" are required to explain 
the observed abundance of the organic refractory dust component in the 
diffuse interstellar medium. Laboratory grain-mantle analog studies 
have provided good experimental evidence for an efficient return of 
grain-mantle molecules to the gas phase upon warm-up of photolyzed 
mixtures, as well as for the formation of an organic refractory residue 
(Hagen et al. 1979; d'Hendecourt et al. 1982; d'Hendecourt 1984; 


ing the accretion time-scale for gas-phase species (eq. [10]). At a 
density of 10^ cm"3 this time-scale is about 10^ yr, which is less than 
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d'Hendecourt et al. 1986; van de Bult 1986). Simple theoretical argu- 
ments, such as those in §8 and below, suggest that UV photolysis and 
transient heating of grains can indeed be of importance. 

In this section, the physical processes that may return grain- 
mantle molecules back to the gas phase are discussed. Some observa- 
tional signatures of grain chemistry will also be indicated. It should 
be emphasized that there is abundant evidence for the importance of gas- 
phase reactions in the formation of simple gas-phase molecules. In 
particular, the molecular ions, HCO + and DC0 + , provide key support for 
ion-molecule gas-phase reactions. 

9.1 Ejection of Grain-Mantle Molecules 

If the accretion of gas on grains is offset by ejection of the accreted 
material, then grains may play an important role in the composition of 
the gas phase. Model calculations of the gas-phase composition of 
molecular clouds in which grain-surface reactions were used alone or in 
conjunction with gas-phase reactions have been performed by several 
groups (Allen and Robinson 1977; Viala et al. 1979; d'Hendecourt et al. 
1985). The contribution of grains to the gas-phase composition in all 
of these models hinges to a large extent on the identification of ejec- 
tion mechanisms for molecules formed on the grain surface. Among the 
ejection mechanisms suggested are thermal evaporation (including single- 
photon heating events for very small grains), ejection upon reaction, 
photodesorption and release of stored chemical energy (cf. Watson 1976; 
Helens and Hagen 1982; d'Hendecourt et al. 1982). 

The first two processes can be of importance for H 2 because of its 
low binding energy and because of its small mass, which makes energy 
transfer to the surface relatively inefficient (cf. §4.6). For other 
molecules, the efficiency of these ejection processes is expected to be 
low. 

Photodesorption is probably an efficient ejection mechanism for 
small molecules on grain surfaces in the diffuse interstellar medium or 
at the edges of molecular clouds. The absorption of a UV photon by a 
molecule on the surface will electronically excite this molecule. 
Photodesorption may then occur because of the electronic repulsion 
between the excited admolecule and neighboring surface molecules. 
Essentially, the larger orbital distribution of the electronically 
excited state coupled with the asymmetric force field of the surface 
will give the nucleus of the admolecule momentum perpendicular to the 
surface, and ejection may occur (Nishi et al. 1984). For small mole- 
cules the efficiency of this process can be fairly large (Greenberg 
1973). Larger molecules will, however, have a relatively low photo- 
desorption efficiency because many more energy-decay channels are 
available and because they may have a larger binding energy. The 
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photodesorption yield will also be low when the excitation energy can be 
rapidly transferred to other molecules in the particle (e.g., near- 
resonant transfer with a time-scale of the order of 10” '3-10”'* sec). 
Thus, for an H2O molecule on a clean substrate the photodesorption yield 
is much higher than for H 2 0 adsorbed on an H 2 0 surface (Greenberg 
1973). Using values typical for small molecules and assuming the aver- 
age interstellar UV radiation field, photodesorption can keep a grain 
surface clean for A y less than about 3 mag (Helens and Hagen 1982). 

The most important mechanism for ejection of molecules in a grain 
mantle into the gas phase inside a dense molecular cloud is probably the 
release of stored chemical energy. At low temperatures, radicals pro- 
duced by UV photolysis can be stored in a grain mantle up to a concen- 
tration of about \% (cf. §8). Upon warm-up, these radicals can diffuse 
and react upon encounter. Laboratory studies on warm-up of photolyzed 
mixtures show that upon heating to about 27 K, those samples with a poor 
heat contact with the substrate will explosively eject most of the : 
mixture (d'Hendecourt et al. 1982; d'Hendecourt 1984). The critical 
temperature of 27 K seems to be a general property of irradiated molecu- 
lar ices, for it is independent of the mixture. Apparently, at this 
temperature rapid diffusion of most of the trapped radicals occurs 
simultaneously. 

It 3hould be emphasized that these explosions are not a result of a 
chain reaction (e.g., the reaction heat of a recombining radical driving 
the diffusion and reaction of other radicals). The deposition of 5 eV 
of reaction heat in a volume corresponding to 200 molecules (i.e., a 
radical concentration of 1?) raises the temperature to about 170 K (see 
below). The temperature relaxation time T, owing to heat diffusion, is 
given by (Anderson 1 98 1 ) as 

t t = 4pCr^/(A) , (18) 

where 0 is the density; C the specific heat 

(»2 x 10^ T 1 '^ erg cm"^ K; Leger et al. 1985); r Q the radius of the 
initial deposition region (-12 A); and K the thermal conductivity 
(a3 x 10 4 erg sec” cm” K ; d'Hendecourt 1984). Inserting these 
values yields a temperature relaxation time-scale of 3 x 10 -1 ^ sec. 
Within that time-scale, however, a radical can migrate at best only a 
few sites and, given the low concentration of radicals, it will not have 
found a reaction partner. Thus, the heat liberated by the reaction does 
not promote reaction of other stored radicals. 

It is probable that the existence of a unique critical temperature 
is a general property of amorphous molecular matrices. At about 25 K 
the heat conductivity of an amorphous material changes abruptly in 
character (Anderson 1981). B^low th^t temperature the energy is carried 
by low-frequency phonons (<10 sec” ) which have a long mean free path 
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( a 1 0^ A). These phonons interact only with tunneling sites where the 
amorphous solid can tunnel from one amorphous state to another (cf. 

§2). Above that temperature the heat is carried by high-frequency 
phonons (sio sec" ) with a mean free path of about 10 A. It is these 
latter, localized phonons that assist in diffusion processes. 

We will now estimate the amount of grain-mantle material evaporated 
by one cosmic ray event or a grain-grain collision. We will assume that 
enough energy is deposited to.raise the grain temperature to the criti- 
cal temperature (AE 3 5 x 10” erg for a 1000 A grain at 27 K) . For 
simplicity it is assumed in this calculation that the grain consists 
mainly of H 2 0 with a trace ( 1% ; cf. §8) of radicals that will start to 
diffuse and react at this temperature. The reaction energy liberated 
raises then the temperature enough to cause appreciable evaporation 
of the mixture. With a concentration of radicals, assuming a 
reaction energy of 5 eV and a specific heat of H 2 0 ice of about 

2 x 10^ T 1 *^ erg cm"^ (T > 50 K); Leger et al. 1985), this yields a 
temperature of about 170 K. This is well above the evaporation tempera- 
ture of H 2 0 ice in the interstellar medium (90 K; Nakagawa 1980), and 
evaporation of the mantle will dominate over other heat-loss mechanisms, 
such as radiation. The amount of material evaporated can then be esti- 
mated by equating the reaction heat with the heat loss through evapora- 
tion. The energy remaining in the particle at the evaporation tempera- 
ture should be subtracted since this energy will be carried away by 
radiation. For the assumptions stated above, this leads to 0.019 eV per 
H 2 0 molecule available for evaporation, which should be compared to the 
binding energy of 0.05 eV. Thus, about 40* of the H 2 0 molecules will 
evaporate before the temperature will drop below the evaporation temper- 
ature. Note that this is less than estimated by Leger et al. (1985), 
because they assumed that grain mantles consist mainly of CO with a 
binding energy of only about 0.09 eV. Cosmic ray heating seems to 
dominate the transient grain heating process with a time-scale of about 

3 x 10^ yr for a grain size of 1000 A. Grain-grain collisions are only 
important at high densities (n Q > 10° cm"'*). With the efficiency calcu- 
lated above, the ejection time-scale for a molecule is thus about 

105 yr. 

9.2 Signatures of Grain Chemistry 

Molecules produced by grain-surface chemistry are typically neutral, 
saturated, small, and highly hydrogenated. Infrared absorption studies 
of protostars show that H 2 0 and alcohols (presumably CH^OH) dominate the 
composition of interstellar grain mantles; NH^ is also present as a 
trace (Tielens et al. 1986a). These molecules are the grain-mantle 
reaction products of the abundant gas-phase species 0, 0 2 , CO, and N. 

The measured abundances of H 2 0 and CH^OH in grain mantles ranges from 
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4 x 10“® to 2 x 10" 1 * per hydrogen atom, which is much larger than their 
gas-phase abundances. Obviously, these molecules, which are easily made 
on grain surfaces, would be prime candidates for studies of the influ- 
ence of grains on the gas-phase composition. In this respect, it should 
be noted that species produced by gas-phase chemistry are often highly 
unsaturated, radicals, or ions. 

In general, the relative contribution of grains to the composition 
of the gas phase is expected to increase with density. Gas-phase chem- 
istry is driven by ion-molecule reactions. The degree of ionization 
decreases, however, with increasing density, essentially because recom- 
bination increases faster with increasing density that ionization 
(Oppenheimer and Dalgarno 1974; Prasad and Huntress 1980). Thus, the 
rate of ion-molecule reactions is expected to increase slower with 
density than the rate of grain-surface reactions. The rate of produc- 
tion of species on grain surfaces is limited by the accretion rate of 
the appropriate reactants onto grains. At molecular cloud densities, 
this means that only for those molecules that are easily made on a grain 
surface but are slowly produced in the gas phase, grains will contribute 
to the gas-phase abundance. 

Thus, one can test the importance of grain chemistry for the compo- 
sition of the gas phase by searching for gas-phase molecules in molecu- 
lar clouds that are difficult to form by ion-molecule reaction 
schemes. Obvious candidates are molecular species containing metal 
atoms (e.g., NaH or NaOH). These species are not expected to form 
molecular species through gas-phase reactions (cf. Oppenheimer and 
Dalgarno 1974; Smith et al. 1983). However, these species will react on 
grain surfaces. 

One might expect, naively, that the overabundance of atomic hydro- 
gen will convert these species into hydrides, as observations show has 
happened to oxygen, nitrogen, and carbon monoxide; however, the grain- 
surface chemistry of such species is somewhat more complicated. In 
particular, hydrogen abstraction reactions can be of importance • (Barlow 
1978). For example, the reaction HaH H ♦ Na ♦ H 2 is exothermic by 
about 2.5 eV. Because a molecular bond has to be broken in order for 
this reaction to proceed, there will be an activation barrier that has 
to be overcome. Typically, such an activation barrier is about 5% of 
the bond energy (*1000 K), which is of little consequence on grain 
surfaces. Thus, in effect, Na like N 2 H^, H 2 C0, and other molecules will 
act as chemical bonding sites for H 2 formation (cf. §5). Hydrogen 
addition and abstraction reactions will continue to form H 2 with Na, 
until it reacts with another migrating species such as 0, which effec- 
tively poisons the site for H 2 formation. Since the energy barrier i 3 
expected to be low and because atomic oxygen is the next most abundant 
migrating radical on a grain surface, most sodium atoms are expected to 
wind up in NaOH molecules. Hydrogen abstraction from these molecules is 
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presumably inhibited by a high actuation energy. It is not yet possible 
to search for the vibrational transitions of molecules with such low 
abundances in interstellar grain mantles, but it is possible to search 
for the rotational transitions of the gas-phase species. 

Nitrogen and sulfur bearing, gas-phase molecules (e.g., NH^ and 
H 2 S) seem difficult to form by ion-molecule reactions (Prasad and 
Huntress 1980, 1982). This would also make them prime candidates for 
studies of the importance of grains for the gas-phase composition of 
molecular clouds. Both NH^ and H 2 S have tentatively been identified in 
grain mantles (Knacke et al. 1982; Geballe et al. 1985). However, the 
gas-phase chemistry of these compounds is still open to some debate (cf. 
Graedel et al. 1982), which somewhat limits their use as grain-chemistry 
indicators. Complex molecules seem also to be difficult to make by gas- 
phase reactions. Grain-surface reactions will not be much help because 
they generally lead to fairly simple, completely hydrogenated mole- 
cules. Ultraviolet photolysis of grain mantles will, however, convert 
such molecules into more complex ones which may eventually be injected 
into the gas phase. The importance of this process is difficult to 
assess. The first step seems to be to identify such molecules in inter- 
stellar grain mantles. 

As discussed before, grain-surface chemistry can lead to large 
deuteration effects (cf. §7 and Fig. 14). Essentially, this reflects 
the high gas-phase ratio of atomic D to H. This ratio is high because 
of the rapid destruction of HD through reactions with H^ (Tielens 1983b; 
cf. §7). Generally, high density and low gas-phase temperatures lead to 
high abundances of deuterated molecules with respect to their hydroge- 
nated counterparts. The temperature dependence arises from the tempera- 
ture sensitivity of the gas-phase H 2 D + plus H 2 reactions. Prime candi- 
dates for such studies are HDO, NH 2 D, DNO, HDCO, and possibly CHjDOH. 
Studies of the deuteration of ammonia may be particularly revealing, 
since the ratios NH 3-i D i/ NH 2-i^i>1 (i = 0, 1,2) should all be equal to 
the gas-phase n(H)/n(D) ratio, which at low temperatures is Just 
A(D)n , where A(D) is the elemental abundance of deuterium. 

Finally, the ortho/para ratio can serve as a grain-chemistry indi- 
cator. Some molecules exist in two distinct states, depending on the 
spin of the nucleus: the para state when the spins are anti-parallel 

and the ortho state when they are parallel. Because of the symmetry of 
the wave function, the para states are associated with the even rota- 
tional levels and the ortho states with the odd levels. Rotational 
transitions between ortho and para states require a simultaneous transi- 
tion of the nuclear spins of the molecule. The radiative lifetime for 
decay from the ortho to the para state is therefore long. The energy 
difference between these two states can be quite large compared to 
interstellar temperatures (i.e., for H 2 about 170 K and for H 2 0 about 
55 K), and small variations in temperature can thus have large effects. 
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The initial ortho-para ratio of molecules such as H 2 0, NH^, H 2 CO, 
and CHjj, newly formed on a grain surface, is expected to be equal to the 
high-temperature limit (i.e., the ratio of the statistical weights of 
the two nuclear spin states). This is due to the large reaction heat 
available upon formation and the considerable rearrangement of the 
hydrogen nuclei involved in most reactions. Ortho-para transformation 
will occur on a time-scale of a day because of the interaction of the 
nuclear spins of neighboring molecules. An even faster ortho-para 
transformation time-scale may result from the presence of paramagnetic 
impurities or lattice defects (Tielens 1986b). The nuclear spin transi- 
tion is due to the inhomogeneous magnetic dipole field generated by the 
nuclear spin of a neighboring molecule or the unpaired electron of a 
paramagnetic impurity or lattice defect. This couples to the proton 
magnetic moments in such a way as to cause a transition between the 
ortho and para states (Wigner 1933). Since the residence time of mole- 
cules produced on grain surfaces is expected to be about 10® yr, 
their ortho-para ratio will reflect the dust temperature (i.e., 
n(o)/n(p) = g exp(-E Q p/T d ) , where g is the ratio of the statistical 
weights of the nuclear spin states, E Q p the ortho-para energy differ- 
ence, and T d the dust temperature). Because ejection into the gas 
phase is a rapid process, it will not influence this ratio. 

Tt should be emphasized that once the molecules are in the gas 
phase, ion-molecule reactions may influence the ortho-para ratio. 
Probably the most important reaction in this respect is H exchange with 
H*. Assuming a Langevin reaction rate coefficient (10*^ cm^ sec” 1 ) and 
a relative abundance of H* of 10*® (Smith and Adams 1984; Lepp et al. 
1986) yields a time-scale ror this process of 3 * 109/n Q yr- This ia 
comparable to the time-scale for accretion on grains (10^/n Q yr), which 
is probably the limiting step in the grain-surface reactions. 

From this discussion, it is clear that in a region where the dust 
and gas temperature are sufficiently different relative to the ortho- 
para energy difference, the ortho-para ratio can be used as a discrimi- 
nant for grain-gas interaction (Tielens 1986b). It should be emphasized 
that the ortho-para transformation process discussed above is not a 
chemical mechanism; strictly speaking, the ortho-para ratio will there- 
fore not prove the importance of grain chemistry for the composition of 
the gas phase. However, the importance of grain chemistry for the 
composition of grain mantles is already well attested to by IR 
spectroscopy . 

The best objects in which to study the ortho-para ratio are cold, 
dense, dark clouds such as Taurus. One such study (for H 2 C0) has been 
performed for two dark clouds (and the warm gas in Orion; Kahane et al. 
1984). Since the energy difference between the ortho and para states of 
H 2 C0 is small (15 K) and the lines suffer from beam dilution, the inter- 
pretation is not straightforward; nevertheless it seems that in dark 
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clouds thermal izat ion on grains dominates, implying an efficient gas- 
grain exchange mechanism. 
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